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Preface 


This  thesis  was  intended  to  determine  whether  airblast 
energy  from  an  atmospheric  nuclear  burst  could  effectively 
couple  into  the  ground,  generate  surface  seismic  signals, 
and  thus  independently  cause  or  contribute  to  the  formation 
of  an  airshock  precursor.  The  major  thrust  of  the  work  was 
to  investigate  the  generation,  description,  and  timing  of 
these  surface  seismic  signals. 

Appreciation  is  extended  to  Major  George  Nickel  of 
the  AFIT  Physics  Department  for  his  recognition  of  an  effort 
worthy  of  research  and  for  his  earnest  attempts  to  teach  free 
thinking  and  the  independent  analytical  application  of  the 
first  principles  of  physics. 


Richard  N.  Price 
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Abstract 


The  coupling  into  the  ground  of  airblast  energy  from 
an  atmospheric  nuclear  burst  is  postulated  as  a mechanism 
which  may  contribute  to  if  not  independently  cause  the 
observed  airshock  precursor.  A computer  model  to  test  the 
hypothesis  is  constructed  by  assuming  an  elastic  ground 
medium,  applying  finite  difference  techniques  to  the  equa- 
tions of  motion,  and  using  the  space-  and  time-varying 
overpressure  from  the  nuclear  burst  to  induce  the  seismic 
motions  within  the  ground. 

The  surface  velocities  resulting  from  simulation  of  a 
28  kiloton  atmospheric  burst  at  500  feet  height  of  burst 
yielded  a dust  layer  ballistically  reaching  only  0.6k  cm  at 
its  highest  point  for  the  stiff  one-layer  ground  medium, 

0.096  cm  for  the  softer  one- layer  medium,  and  a negligible 
height  for  the  more  realistic  four- layer  Frenchman  Flats 
medium.  Thus,  the  airblast-induced  precursor  as  postulated 
(ballistic  rise  only)  fails  to  re-create  the  2-3  meter  high 
dust  layers  observed  in  experimental  atmospheric  nuclear 
testing.  However,  the  motions  are  felt  to  be  significant 
enough  to  be  included  in  any  attempt  to  model  from  first 
principles  the  precursor  and  the  up-sweep  of  dust  behind  the 


shock  front 


I. 


Introduction 


Background  and  Phemonenology 

During  early  atmospheric  nuclear  weapons  testing 
the  effect  known  as  the  nuclear  airhlast  precursor  was 
observed.  The  precursor  results  from  the  formation  of 
a heated  layer  of  air  immediately  above  the  ground  surface. 
Because  the  sound  speed  is  higher  in  this  heated  air  layer, 
the  shock  front  from  the  nuclear  burst  is  able  to  propagate 
faster  in  this  heated  region  than  in  the  higher,  cooler 
regions  of  air.  Thus,  a "toe"  forms  on  the  leading  edge  of 
the  shock  front  at  its  intersection  with  the  ground.  This 
phenomenon  is  shown  pictorially  in  Figure  1 . 

A nuclear  burst  which  generates  an  airblast  precursor 
represents  a departure  from  the  shock  front  shape  and  shock 
properties  one  would  expect  when  a nuclear  burst  is  detonated 
over  an  ideally  reflecting,  non-interacting  surface.  It  is 
precisely  the  interaction  of  a real  surface  which  gives  rise 
to  the  non-ideal  (precursed)  properties  of  the  airblast: 
two  overpressure  peaks  instead  of  one;  peak  overpressure 
reduced  up  to  50%>  over  ideal;  and  dynamic  pressure  increased 
up  to  100%>  over  ideal.  The  properties  of  a nuclear  blast 
precursor  therefore  become  of  considerable  interest  when 
the  response  of  surface  and  near -surface  systems  to  the 
static  loading  of  real-world  overpressure  waveforms  and  to 
the  enhanced  drag  loading  from  the  dynamic  pressure  must  be 
known  with  reasonable  accuracy. 
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Altitude *-  Altitude 


The  Air  Force  Weapons  Laboratory  (AFWL)  began  to 
study  the  precursors  observed  in  actual  atmospheric  nuclear 
tests  (Table  I)  in  an  attempt  to 

1.  characterize  the  development  of  the 
heated  air  layer  which  is  the  heart  of 
the  precursor; 

2.  model  the  development  of  this  thermal 
region;  and 

3.  extend  this  understanding  of  experimental 
observations  for  the  relatively  low  yield 
experimental  nuclear  blast  precursors 

(1  to  30  kilctons)  to  the  megaton  yields 
which  the  defense  planner  must  expect 
in  today's  threat  environment. 

Utilizing  considerable  internal  and  contractor 
expertise  from  1973  to  the  present.  AFWL  has  principally 
concluded  that  radiation  (primarily  X-rays)  from  the  fireball 
is  the  mechanism  by  which  the  heated  air  layer  ( and  thus 
the  precursor)  is  formed.  Specifically,  if  the  thermal 
radiation  on  the  ground  is  sufficient  to  heat  the  soil  above 
a threshold  level  (which  will  depend  on  soil  type  and  soil 
conditions),  ground  moisture  and  hydrated  water  in  the  soil 
will  be  suddenly  and  explosively  released.  Such  a violent 
release  will  "popcorn"  the  soil  into  the  air.  Once  airborne, 
the  dust  heats  the  surrounding  air  by  conduction  and  convec- 
tion while  continuing  to  absorb  thermal,  radiation  from  the 
fireball.  (Ganong,  1978) 

Nonetheless,  the  coupling  of  airblast  energy  to  the 
ground  has  been  postulated  as  another  causal  mechanism  which 
may  contribute  to  if  not  independently  cause  the  precursor. 
This  mechanism  has  not  been  previously  studied  to  determine 


3 


(J 


I 


Table  I.  Events  Studied  for  Precursor 


Shot 

Operation 

(Number) 

Yield 

(Kt) 

HOB 

(ft) 

Scaled 

non, 

(ft  ■/Kt  1/3) 

Area 
al  NTS 

Type 

POST 

Teapot.  (11) 

1.45 

300 

265 

9 

Tower 

MOTH 

Teapot  (2) 

2.40 

300 

224 

3 

Tower 

HORNET 

Teapot  (5) 

3.6 

300 

196 

3 

Tower 

TESLA 

Teapot  (3) 

6.8 

300 

158 

9 

Tower 

FOX 

Tumbler-Snappei  (G) 

11.5 

300 

133 

4 

Tov/er 

EASY 

Tumbler-Snapper  (5) 

12.  5 

300 

129 

1 

Tower 

HO'V 

Tumbler-Snapper  (8) 

13.9 

300 

125 

3 

Tower 

NANCY 

Upshot-Knothole  (2) 

24 

104 

4 

Tower 

BADGER 

Upshot-Knothole  (G) 

25 

1C3 

2 

Tower 

HARRY 

Upshot-Knothole  (8) 

32.3 

300 

94 

2 

Tower 

SIMON 

Upshot -Knothole  (7) 

45 

84 

1 

Tower 

MET 

Teapot  (12) 

22.5 

400 

142 

FF 

Tower 

BEE 

Teapot  (6) 

8 

500 

250 

7 

Tower 

MORGAN 

Plumbbob  (3b) 

8 

250 

9 

Balloon 

V/ILSON 

Plumbbob  (5) 

10.3 

230 

9 

Balloon 

KEPLER 

Plumbbob  (11) 

10.3 

230 

4 

Tower 

BOLTZMANN 

Plumbbob  (2) 

11.5 

1 1 

222 

7 

Tower 

APPLE  I 

Teapot  (8) 

14.2 

500  1 

206 

4 

Tower 

SHASTA 

Plumbbob  (16) 

10.5 

•coo 

196 

2 

Tower 

WHITNEY 

Plumbbob  (28) 

18.5 

500 

189 

2 

Tower 

ZUCCHINI 

Teapot  (14) 

28 

500 

165 

7 

Tower 

APPLE  II 

Teapot  (13) 

28.5 

500 

164 

1 

Tower 

TURK 

Teapot  (4) 

44 

500 

142 

2 

Tower 

GRABLE 

Upshot-Knothole  (1G) 

15 

524 

212 

FF 

Gun 

PRISCILLA 

Plumbbob  (G) 

3G.G 

700 

211 

FF 

Balloon 

SMOKY 

Plumbbob  (20) 

44 

700 

198 

8 

Tower 

WASP  PRIME 

Teapot  (9) 

3.2 

■9 

502 

7 

Air 

LA  PLACE 

Plumbbob  (24) 

1.22 

E3 

702 

7 

Balloon 

V/ASP 

Teapot  (1) 

1.2 

7G2 

717 

7 

Air 

ABLE 

Tumbler-Snapper  (1) 

1 

793 

793 

FF 

Air 

DOG 

Tumbler-Snapper  (4) 

13.5 

1040 

393 

7 

Air 

CHARLIE 

Buster- Jangle  (3) 

14 

1132 

470 

7 

Air 

CLIMAX 

Upshot- Knothole  (1 1) 

GO 

1334 

341 

7 

Air 

DOPPLER 

Plumbbob  (17) 

19.7 

1500 

631 

7 

Balloon 

NEWTON 

Plumbbob  (2G) 

11.  b 

1500 

u59 

7 

Balloon 

! STOKES 

Plumbbob  ( 14) 

19 

1500 

562 

7 

Ba'lnon 

Precursor1 


Yes 


Yes 

Yes 


No 

? 

Yes 

Yes 

Yes 

""Ygs"" 


Yes 


Yes 

Yes 

Yes 
? 

? 

? 


(Liner,  1975s  22) 


its  relative  contribution  to  the  heated  air  layer  which 
precedes  and  determines  the  precursor  formation. 

The  airblast-generated  precursor  contribution  as 
postulated  would  occur  in  the  following  sequence  (Figure  2): 
the  airblast  strikes  the  ground,  coupling  energy  into  the 
ground  and  generating  seismic-like  displacements;  the  dis- 
placement waveforms  propagate  into  the  ground  along  the 
soil  surface  with  velocities  characteristic  of  the  seismic 
velocities  of  the  ground  medium.  As  the  radius  of  the 
airshock  increases,  the  peak  overpressure  at  the  shock  front 
decreases  and  the  airshock  subsequently  slows.  When  the 
seismic  surface  wave  velocity  exceeds  the  velocity  of  the 
advancing  airblast,  these  seismic  surface  disturbances  will 
outrun  the  airshock  and  are  for  the  first  time  able  to 
contribute  to  tne  thermal  layer  formed  ahead  of  the  shock. 
For  this  contribution  to  be  of  concern,  the  surface 
disturbances  must  be  sufficiently  large  to  impart  a sizable 
vertical  velocity  to  the  dust/soil  particles  lying  loosely 
on  the  ground,  causing  these  loose  particles  to  rise 
ballistically  above  the  ground.  One  additional  requirement 
is  that  sufficient  energy  must  remain  in  the  fireball  to 
heat  the  rising  dust  particles.  The  heated  dust  can  then 
heat  the  surrounding  air  by  conduction  and  convection  as 
in  the  case  of  the  thermally-induced  soil  blowoff.  Early 
heacing  from  the  fireball  will  practically  assure  the 
presence  of  loose,  dry  soil  at  the  ground  surface.  As  in 
the  case  of  thermally- induced  soil  blowoff,  the  airblast- 


Figure  2 


! 


* 


a.  Airblast  strikes  ground  at 


Tq;  Fireball  heats  ground 
surface. 


* 


b.  Airblast  expands  along  sur- 
face; Seismic  waves  radiate 
into  ground  and  along  surface 
but  move  slower  than  front  of 
airshock. 


c.  Seismic  waves  are  now  able  to 
outrun  airblast,  imparting  an 
upward  velocity  to  dust  parti- 


induced  contribution  will  also  be  yield,  height  of  burst, 
and  soil  condition  dependent. 


Purpose 

The  documentation  which  follows  details  the  approach 
that  was  undertaken  to  model  any  airblast-generated  contri- 
bution to  the  heated  air  layer  which  precedes  the  precursor. 
As  such,  it  presents  computations  of  seismic  surface  dis- 
placements, seismic  surface  velocities,  and  anticipated 
airborne  dust  layers.  This  research  is  intended  to  be 
used  to  determine  whether  further,  more  exact  modeling  is 
warranted  on  the  role  of  airblast  energy  couples  into 
seismic  waves  as  a contributor  to  the  nuclear  precursor. 

Scope 

The  computational  model  to  simulate  ground  motions 
developed  as  a product  of  this  thesis  research  is  applicable 
to  any  problem  which  involves  pressure  or  stress  loading 
normal  to  the  surface  of  an  elastic  half-space  for  which 
order  of  magnitude  answers  are  desirable.  However,  the 
problems  for  which  ground  motion  results  are  presented  are 
limited  to  atmospheric  nuclear  test  events.  This  work  does 
not  investigate  thermal  energy  transport  or  hydrodynamics 
within  the  rising  dust  layer. 


■ 
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Order  of  Presentation 


The  order  of  presentation  will  be  as  follows i 
Chapter  II  describes  the  computational  model  (preliminary- 
analysis,  equations  of  motion,  finite  difference  equations, 
surface  pressure,  treatment  of  boundary  displacements,  and 
method  of  solution) ; Chapter  III  discusses  output  calculated 
from  the  model  (stability,  transmitting  boundary,  and 
simulation  of  an  atmospheric  nuclear  burst);  and  Chapter  IV 
concludes  whether  the  seismic  surface  waves  generated  by 
the  airblast  make  significant  contributions  to  the  precursor 
and  recommends  improvements  or  new  approaches  for  the  model. 


II. 


Computational  Model 


Preliminary  Analysis 

The  heart  of  the  airblast  precursor  is  the  formation  of 
a heated  air  layer  immediately  above  the  ground.  The  postulated 
means  through  v/hich  seismic  signals  can  contribute  to  this 
thermal  layer  is  by  the  injection  of  dust/soil  particles  into 
the  air,  these  particles  subsequently  absorbing  fireball 
radiation. 

The  height  h to  which  a dust  particle  would  rise 
ballistically  above  the  surface  (ignoring  hydrodynamics  and 
drag)  is  independent  of  its  mass  and  is  given  by 


where  v = the  particle's  vertical  velocity,  and 
g = the  gravitational  acceleration. 

Clearly  then,  the  driving  ^oree  involved  in  tasting 
this  postulated  causal  mechanism  is  to  calculate  an  order  of 
magnitude  value  for  the  vertical  velocity  of  the  dust/soil 
particles,  or  equivalently,  the  peak  vertical  velocity  of 
the  ground  surface.  The  single  greatest  assumption  made  is 
that  satisfactory  values  can  be  obtained  through  a treatment 
of  the  soil  as  an  elastic  half-space. 

This  assumption  is  put  into  perspective  through 
consideration  of  research  conducted  by  the  U.  S.  Army  Corps 
of  Engineers.  Since  1971  computations  and  data  analyses  have 
simulated  and  reconstructed  seismic  motions  resulting  from 


high  explosive  and  nuclear  test  events.  This  work  indicates 
that  while  soils  are  not  truly  elastic  media,  two-dimensional 
elastic  wave  propagation  calculations  can  be  performed  with 
comparative  ease  and  require  only  limited  soil  property 
data  relative  to  inelastic,  non-linear  calculations.  Also,  it 
was  felt  that  elastic  calculations  were  less  subject  to 
numerical  errors.  (Hadala,  1973*  297,  383) 

In  addition,  one  particular  set  of  computations  compiled 
by  the  Army  Corps  of  Engineers  revealed  that  for  peak  over- 
pressures in  the  100  to  50  psi  range,  the  elastic  soil  model 
generally  yielded  maximum  velocities  less  than  those  given  by 
the  inelastic  model.  However,  the  reverse  was  true  in  the  50 
to  1 psi  range.  In  other  related  calculations  a similar  but 
different  overpressure  crossover  point  was  also  found.  Generally 
then,  these  results  imply  that  the  elastic  model  subsequently 
used  to  approximate  the  vertical  velocities  of  the  dust 
particles  will  underestimate  these  velocities  if  outrunning  of 
the  airshock  by  the  seismic  signal  occurs  at  the  high  end  of 
the  peak  overpressure  region  and  will  overestimate  them  if 
outrunning  occurs  in  the  low  overpressure  region.  (Hadala,  1973* 


295) 


After  first  assuming  an  elastic  medium,  an  estimate 


through  one-dimensional  analysis  can  be  made  of  the  maximum 
surface  vertical  velocity  in  a uniform  media.  This  is  derived 


in  Appendix  A and  is  given  by 


d Sz 


max  at  z=0 


ground  ^c 


(cm/sec) , 


where  S = displacement  in  the  z-direction, 
z 

P = peak  overpressure  on  the  surface  at  the 
0 point  of  interest, 

ground  = density  of  the  ground,  and 

V = compressional  seismic  velocity  of  the 
c ground. 

The  peak  overpressure  Pq  along  the  surface  is  a function 
of  yield  and  height  of  burst;  and  the  compressional  seismic 
velocity  Vc  is  a function  of  the  seismic  constants  (Lame/ 
constants  \ and  jU  ) and  density  of  the  ground.  Therefore, 
the  maximum  vertical  velocity  with  which  dust  can  be  injected 
into  the  air  will  also  be  dependent  upon  the  yield,  height  of 
burst,  and  seismic  properties  of  the  ground. 

Another  important  consideration  - the  ground  range  at 
which  the  seismic  signal  will  begin  to  outrun  the  airshock  - 
will  also  have  the  same  dependency  as  above.  In  order  to 
obtain  a preliminary  estimate  of  where  this  "outrunning"  will 
occur,  it  is  fi^st  assumed  that  there  is  an  average  velocity 
Vave  with  which  the  sperical  shock  expands.  With  this 
assumption  it  is  shown  in  Appendix  G that  the  apparent  velocity 
with  which  the  airshock  advances  along  the  ground  surface  can 
be  expressed  as  f i 


(cm/sec)  (lb) 


2 2 
n + r 


where  H = height  of  burst,  and 
r = ground  range. 


When  this  velocity  V 


slows  to  the  seismic  velocity 


app 

V , the  seismic  signal  thus  created  at  the  intersection  of  the 

v 

shock  front  and  ground  will  outrun  the  airshock. 

However,  the  surface  seismic  disturbances  having  the 

largest  magnitude  are  known  to  be  Rayleigh  waves  which  travel 

more  slowly  than  V . So,  solving  Equation  lb  for  r when  V = 

c ^ app 

VRayleigh’  y^-elds  that  ground  range  at  wnich  outrunning  of 
seismic  Rayleigh  waves  will  occur.  This  is  given  by 


H 


V 

Rayleigh 

V 

L ave 


- 1 


(cm  or  ft) . 


Consider  the  case  of  a 30  kiloton  burst  detonated  at 
500  feet,  also  where  V = 2.4  km/sec  (eight  times  the  sound 
speed  in  sea-level  air)  and  = 2.5  km/sec.  Then, 

Rayleigh  outrunning  begins  to  occur  at  a ground  range  of 
r = O.52  km  = 1714  ft.  From  blast  data  curves  it  is  found 
that  at  that  ground  range  the  peak  overpressure  will  be 
PQ  ==  40  psi  and  that  the  time  of  arrival  of  the  airshock  will 
be  TARR  ~ 0.15  sec.  Thus  the  maximum  vertical  velocity 
(Equation  la)  is  2.6  cm/sec  and  50  percent  of  the  fireball 
radiation  has  yet  to  be  emitted.  (Glasstone,  1977:  111,  121, 

309-310) 

While  Equations  la  and  lb  were  derived  for  a uniform 
elastic  medium,  peak  velocities  occurring  in  a layered  elastic 
medium  are  also  of  interest  because  ground  motion  data  and 


12 


soil  testing  have  revealed  that  testing  grounds  such  as 
Frenchman  Flats  at  the  Nevada  Test  Site  consist  of  horizontal 
layers  which  vary  with  depth  in  density  and  characteristic 
seismic  velocities  (Figure  3)*  Because  elastic  property 
data  for  the  medium  at  Frenchman  Fiats  are  known  and  atmospheric 
nuclear  tests  displaying  an  airblast  precursor  were  conducted 
there,  the  layered  medium  of  this  site  will  provide  a "real- 
world"  problem  with  which  to  test  the  postulate  of  the  airblast 
generated  precursor. 

Figure  4 depicts  the  phenomological  effect  of  a two- 

✓ 

layered  ground  medium.  Typically,  because  of  the  effects  of 
weathering,  pressure  from  the  overburden,  and  some  degree  of 
cementation,  deeper  layers  will  be  seismically  "stiffer"  (read 
less  compressive)  and  possess  a higher  characteristic  seismic 
propagation  velocity.  Thus,  once  a disturbance  reaches  a 
stiffer  layer  it  will  propagate  more  quickly  than  a disturbance 
in  an  above  layer.  This  will  ultimately  cause  the  seismic 
signal  to  outrun  the  airblast  at  the  ground  surface  earlier 
than  would  have  occurred  had  the  material  consisted  solely  of 
material  with  seismic  properties  of  the  upper  layer. 

One  other  important  effect  of  layering  in  the  medium 
is  that  the  stress  transmitted  into  the  lower  layer  can  be 
greater  than  the  stress  incident  on  the  interface  from  above. 

The  transmitted  and  reflected  stresses  have  been  derived 
from  one-dimensional  analysis  in  Appendix  B and  are  given  by 
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where 


density  of  layer  j,  and 
seismic  velocity  of  layer  j . 


For  a stiffer  lower  layer,  density  -f2  and  seismic 
velocity  C0  will  be  larger  than  ^ and  of  the  upper  layer. 
Thus  the  transmitted  stress  will  be  greater  than  and  of  the 
same  type  (tensile  or  compressive)  as  the  incident  stress. 
Therefore,  if  the  energy  attenuation  resulting  from  spatial 
expansion  of  the  seismic  signal  is  not  greater  than  the  gain 
in  stress  achieved,  the  stress  transmitted  to  the  second  layer 
which  is  then  later  re-transmitted  back  into  the  upper  layer 
can  be  greater  than  the  original  incident  stress,  thereby 
resulting  in  seismic  displacements  greater  than  in  the  single- 
layered medium. 


Equations  of  Motion 

Armed  with  a knowledge  of  the  seismic  phenomenology 
and  the  parameters  upon  which  it  depends,  the  next  step  is  to 
generate  the  equations  which  govern  the  motion  of  the  ground 
medium  under  airblast  loading.  Cylindrical  geometry  is 
chosen  due  to  the  axial  symmetry  of  the  advancing  airblast 


shock  front.  S and  S are  defined  to  be  the  displacements 
r z 

in  the  radial  and  vertical  directions,  respectively,  of  a 
given  point  within  the  ground  medium  about  its  equilibrium 
point.  The  displacement  Sq  in  the  ©-direction  is  zero  because 
of  the  problem  symmetry. 

The  equations  of  motion  are  given  by 
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The  boundary  conditions  which  apply  are  that  the  stress 
at  the  surface  z=0  in  the  vertical  direction  is  just  the 
overpressure  P(r,t)  which  acts  normal  to  the  surface. 


= A A + 2Ai  = -P(r,t) 


where 


and  that  no  tangential  stress  exists  at  the  surface.  (Kolskys 


195:  55) 


= 0 = 11 


Finite  Difference  Equat j ons 

The  technique  of  finite  differencing  was  chosen  to 
solve  the  equations  of  motion  primarily  due  to  the  ease  and 
speed  with  which  the  equations  could  be  implemented  and  solved 
on  a digital  computer.  The  equations  were  evaluated  at  the 
time-space  mesh  point  i,j,n,  where  the  finite  difference  mesh 
is  described  in  Figure  5 and 

S (r..z.,t  ) = S (iAr,  jAz,  nAt)  = S 

J i jn 

A straight-forward  application  of  central  differences 
yields  for  S^, 


Because  central  differences  were  employed,  the  scheme 

will  be  accurate  to  second  order  in  At,  Ar,  and  Az«  The 

centering  grid  is  given  in  Figure  6. 

Similar  implementation  of  the  finite  differencing  for 

the  S displacement  and  the  boundary  conditions  causes  Equations 
z 

3,  6,  and  7 to  become 

S - 2 S + Sr 

zi>j.n+l  zi,j,n  zi,j,n-l 

A?  = 
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Because  the  displacement  values  at  time  n and  n-1  will 
be  zero  until  airshock  arrival,  the  only  unknowns  in  Equations 
8-11  are  the  time  advanced  (n+l)  values.  Thus,  the  scheme  as 
chosen  is  a second-order  accurate,  explicit  algorithm. 
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Surface  Pressure 
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Finally,  before  implementation  on  the  computer  the 
surface  pressure  input  P(r,t)  is  required.  Because  scaled 
near-proximity  of  the  bursts  to  the  ground  is  necessary  to 
achieve  an  airblast  precursor,  surface  effects  upon  the  over- 
pressure have  to  be  considered.  The  blast  data  curves  in 
Glasstone's,  The  Effects  of  Nuclear  Weapons  were  taken  as 
reference  data  because  of  his  inclusion  of  ground  effects. 

For  the  typical  overpressure  waveform  shown  in  Figure  7, 
the  following  parameters  are  required:  time  of  arrival  TARR, 
peak  overpressure  PEAKP,  duration  of  the  positive  overpressure 
phase  TPLUS,  and  the  rate  of  decay  of  pressure  witli  time. 

General  algorithms  were  found  for  TARR  and  derived  for 
the  rate  of  decay  of  the  pressure  with  time.  The  first  is 
given  as  follows  (Liner,  1975*  1972-3): 


0.54291  Y - 21.185  r Y2^3  + 361.8  r 2Y1//3  + 2383  r 2 

s s s 

TARR  T-77: 5 

Y /J>  + 2.048  r Y1/  i + 2.6872  r 


where 

TARR  is  in  msec, 

Y = weapon  yield  in  kilotons,  and 

r = slant  range  in  kilofeet. 
s 

For  a surface  burst 

TARR  (HOB  = 0)  = TARR  (2Y,  r ). 

s 
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t_  „„  = arrival  time  of  shock  front 

<XL  1 

t*  = duration  of  positive  pressure 
p phase 

(negative  pressure  phase  will  he 
ignored  in  this  model) 
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For  a non-surface  burst,  the  free  air  arrival  time  is 
used  in  the  regular  reflection  region  and  a linear  inter- 


polation between  free  air  and  surface  burst  values  is  made  in 
the  Mach  reflection  region. 

HOB 

TARR  = TARH(Y.rs),  for  range  - li 

HOB 

TARR(Y,rs)  * ground  range  + 


TARR(2Y,r  ) * (l 

s 


HOB 

ground  range 


) , for 


HOB 

ground  range 


The  second  parameter,  decay  of  overpressure  with  time, 
is  given  by 


P(TAU) 

PEAKP 


1 - (1  - TAU)m 


J./m 


where 

TAU 


m 


t - TARR 

TPLUS  ’ and 

1 + .382  (In  PEAKP) 

- .136  (In  PEAKP)2  + .025  (In  PEAKP)3  . 


The  computational  algorithm  for  TARR  was  originally 
derived  by  Brode  and  that  for  P(t)  was  derived  by  applying 
curvefitting  techniques  (Appendix  C)  to  the  Glasstone  data 
for  overpressure  decay.  (Glasstone,  1977*  100) 

At  the  time  of  computer  implementation  no  suitable, 
general  algorithms  had  been  found  for  PEAKP  and  TPLUS.  The 
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technique  chosen,  then,  was  to  fix  The  Glasstone  data  for 
these  parameters  by  linear  segmentation  for  each  specific 
problem  simulated.  This  aspect,  while  speeding  initial 
implementation,  adds  significant  awkardness  when  applying  the 
code  to  various  burst  problems.  A minor  improvement  was 
achieved  by  calculating  the  pressure  as  a function  subprogram. 
For  differing  events  the  unique  overpressure  function  can  be 
validated  independently  and  then  appended  to  the  main  program. 
(Glasstone,  1977s  111-115.  119) 

Treatment  of  Boundary  Displacements 

The  two  stress  components  on  the  surface  (the  normal 
stress  defined  by  the  overpressure  and  the  identically  zero 
tangential  stress)  provide  the  calculational  means  to  derive 
the  displacements  at  the  top  boundary  of  the  finite  difference 
mesh.  These  were  given  analytically  by  Equations  6 and  7 and 
in  finite  difference  form  by  Equations  10  and  11. 

Symmetry  is  the  key  by  which  the  displacements  at  the 
lefthand  mesh  boundary  are  found.  Because  ground  zero  is  the 
symmetry  axis,  displacements  immediately  to  either  side  are 
considered  equal.  By  averaging  these  mirrored  displacements, 
the  displacements  at  the  lefthand  mesh  boundary  or  ground  zero 
are  found  and  are  equal  to  those  just  averaged. 

Finally,  the  righthand  and  bottom  mesh  displacements 
must  be  defined.  The  simplest  choice  is  to  apply  rigid 
boundaries;  that  is,  the  displacements  there  remain  zero.  One 
particularly  annoying  property  of  such  a treatment  is  that 
reflection  of  seismic  signals  occurs  at  these  rigid  boundaries. 


27 


This  trait  becomes  particularly  restrictive  when  attempting 
to  compare  model  output  data  to  published  late-time  seismic 
motions  from  test  data  and  independent  elastic  calculations. 

Interference  of  reflected  signals  can  be  prevented.  Two 
computationally  easy  choices  can  readily  be  implemented: 
increase  the  mesh  size  or  stop  the  computation  before  any 
reflected  signal  can  reach  the  point  of  interest. 

Increasing  the  mesh  size  can  be  done  by  increasing  the 
spatial  mesh  increments,  thereby  losing  accuracy,  or  by 
increasing  the  computer  memory  requirements  of  the  calcula- 
tion. Stopping  the  computation  before  reflection  can  affect 
the  motion  of  the  point  of  interest  can,  as  in  late-time 
motion,  result  in  shutdown  before  significant  motion  has 
occurred  at  the  point  of  interest. 

Another  choice  which  is  more  difficult  computationally 
is  to  develop  a transmitting  boundary;  that  is,  a boundary 
which  acts  as  nearly  as  possible  as  though  a semi-infinite 
region  of  material  exists  beyond  the  mesh  boundary.  If 
successfully  accomplished,  a seismic  signal  can  be  transmitted 
across  the  boundary  with  no  reflection;  and,  the  mesh  could  be 
kept  resonably  compact,  yet  still  yield  acceptably  accurate 
results. 

The  nature  of  the  airblast-induced  precursor  requires 
that  surface  seismic  motions  be  computed  over  a ground  range 
which  is  long  in  comparison  to  the  ground  depth  of  interest. 

As  a result  a relatively  shallow  mesh  grid  is  used.  Because 
reflections  from  the  bottom  of  the  mesh  would  therefore  occur 
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first,  the  bottom  mesh  boundary  was  chosen  for  the  application 
of  the  transmitting  boundary. 

The  requirement  upon  such  a boundary  is  that  it  must 
act  as  a one-way  valve.  This  one-way  action  can  be  accomplished 
by  permitting  a seismic  disturbance  incident  upon  the  boundary 
from  above  to  pass  out  of  the  mesh  and  disappear  into  the 
imaginary  half-space  below  while  simultaneously  preventing 
the  return  of  reflected  signals  into  the  mesh. 

The  means  by  which  this  action  will  be  performed  is  based 
upon  a momentum  flux  argument.  Consider  the  mesh  as  it  exists 
in  Figure  8 near  the  bottom  boundary.  At  the  points  marked 
with  the  symbol  X,  the  displacements  and  are  known.  From 
these  displacements,  calculations  of  velocities,  spatial 
derivatives,  and  stress  components  can  be  made  for  the  point 
P centered  at  i+J  and  JMAX-3/2. 

The  stress  component  0~zz  represents  the  flux  of  z-momentum 
in  the  z-direction,  and  component  q-  , the  flux  of  r-momentum 
in  the  z-direction.  A properly  constructed  transmitting 
boundary  permits  the  net  outward  flow  of  momentum.  However,  it 
prevents  the  net  inward  flow  of  momentum  because  this  results 
in  the  undesirable  increase  of  momentum  within  the  finite 
difference  computational  grid. 

Therefore,  if  the  proper  flow  of  momentum  is  indicated  at 
point  P (Appendix  D presents  the  derivation  of  the  logic  table 
and  equations  for  the  transmitting  boundary) , this  momentum  is 
allowed  to  flow  across  the  boundary  (level  JMAX-1)  by  equating 
the  stress  at  point  P to  that  at  point  R centered  across  the 


29 


boundary  at  i+f-  and  JMAX--|.  If  the  conditions  indicate  improper 
flow  of  momentum,  the  stress  components  at  R are  set  to  zero. 
Next,  displacements  along  level  j= JMAX  are  calculated  from  the 
centering  scheme,  where  the  displacements  at  the  lower  right- 
hand  corner  of  the  mesh  are  zero  and  the  computational  march 
is  executed  along  the  j=JMAX  level  from  i=IMAX-l  to  i=3» 

Of  course,  a trade-off  does  exist  when  using  such  a 
transmitting  boundary.  This  is  that  momentum  transferred  out 
of  the  mesh  is  lost  and  can  play  no  part  in  the  elastic  rebound 
of  the  medium.  For  airblast-induced  seismic  motions  this 
limitation  is  expected  to  occur  far  too  late  to  be  of  concern. 
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Method  of  Solution 

The  explicit,  second-order  finite  difference  scheme 
derived  earlier  involves  two  known  time  levels  (n-1  and  n)  and 
three  known  spatial  levels  in  both  the  radial  ( i— 1 , i,  i+l) 
and  vertical  (j+1,  j , j-l)  directions.  The  task  is  then  to 
calculate  the  unknown,  time-advanced  (n+l)  values  for  all  the 
spatial  mesh  points  of  interest.  Figure  9 presents  as  a flow 
chart  the  method  by  which  the  computational  model  solves  the 
equations  of  motion. 

Data  characterizing  the  burst,  ground  medium,  finite 
difference  mesh  parameters,  and  output  options  are  first  input. 

An  array  describing  the  ground  properties  within  the  mesh  is 
constructed.  Next,  the  overpressures  on  the  ground  surface 
are  calculated  and  used  by  the  boundary  condition  equations. 

These  equations  give  the  displacements  at  the  uppermost 
boundary  of  the  finite  difference  mesh,  and  the  symmetry  condi- 
tion is  used  to  find  the  displacements  at  the  leftmost  boundary. 
Utilising  the  main  difference  equations  for  the  seismic  displace- 
ments and  marching  through  the  spatial  mesh,  the  displacements 
are  found  for  each  mesh  point  except  at  the  righthand  and 
bottom  boundaries. 

No  displacements  are  calculated  at  the  righthand  boundary 
because  a rigid  boundary  condition  is  used.  For  the  bottom, 
either  a rigid  boundary  condition  or  a transmitting  boundary 
condition  based  on  a momentum  transfer  argument  is  applied 
(Appendix  D) . All  of  the  desired  time-advanced  displacements 
throughout  the  mesh  have  not  been  calculated. 
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Figure  9. 


Method  of  Solution 
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Finally,  the  output  format  chosen  by  the  input  options 
is  executed,  and  the  process  is  repeated,  beginning  with  an 
advanced  time  and  updated  surface  pressure. 

This  iterative  process  continues  until  the  desired 
shutdown  time  is  reached.  Typically,  this  shutdown  time  is 
chosen  as  the  time  of  first  arrival  of  seismic  waves  at  the 
righthand  rigid  boundary.  Thus,  unwanted,  artificial  reflections 
are  prevented. 


III.  Code  Calculations 

Stability 

One  of  the  prime  questions  arising  when  using  any 
finite  difference  scheme  is  whether  the  scheme  is  stable. 

The  artificial,  non-physical  inaccuracies  introduced  by  the 
presence  of  instability  are  governed  by  the  choice  of  time 
step  At  and  grid  element  sizes  /\r  and  Az. 

Two  different  stability  criteria  were  initially  speculated 
to  be  applicable  to  this  model.  The  first,  St^ , was  set  forth 
by  the  Army  Corps  of  Engineers  investigators  as  applicable  to 
two-dimensional  finite  difference  calculations  of  ground  shock 
which  outruns  the  airblast.  This  condition  was  presented  as 
the  Courant,  Fredricks,  and  Levy  stability  requirement,  namely, 
(Hadala,  1973s  42) 


0 

c+ 

1 A 

(sec). 

vmax 

where  At  = 

time  step, 

Ax  = 

spatial  mesh  increment  = Ar  = 

A z,  and 

V 

max 

maximum  disturbance  propagation 

velocity. 

The  second  stability  criterion,  St2,  was  based  on  the 
physical  argument  that  no  disburbance  should  be  allowed  to 
propagate  a distance  greater  than  the  mesh  spacing  Ar  or  Az- 
in one  time  step  A't-  Mathematically,  this  is  given  as 

At  *=  (sec) 

~ V 

max 
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In  order  to  determine  which  of  the  two  preceding  stability 
conditions  in  fact  apply  to  this  model,  an  experimental  stability 
analysis  was  conducted.  Figures  10-12  present  the  results  of 
this  stability  analysis.  This  analysis  consisted  of  making 
several  runs  with  the  computer  model  to  calculate  a ground 
shock  problem.  All  runs  are  the  same  except  that  the  time  step 
has  been  changed  in  each  run.  The  time  step  for  each  run  is 
shown  in  relation  to  both  stability  criteria,  St^  and  St^. 

While  some  far-field  oscillations  appear  for  At  = .002  sec 
(which  exceeds  the  stability  criterion  St  set  forth  by  Hadala) , 
the  solution  has  not  become  unstable.  Instead,  this  indicates 
that  the  stability  condition  is  being  approached.  Further 
enlargements  of  the  time  step  show  that  the  less  restrictive, 
physically-derived  condition,  St2,  should  be  applied  as  the 
stability  limitation  to  the  difference  scheme  of  this  model. 
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Figure  10.  Cont. 


Figure  11.  Cont. 


Figure  12.  Stability  Analysis , 2. 5 Msec  Time  Step 
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Transmitting  Boundary 

Figures  13-16  present  a comparison  of  surface  seismic 
motions  computed  by  the  computer  model.  Three  different  but 
related  problems  are  represented.  One  is  a "standard" 
employing  a rigid  bottom  boundary  and  100  by  25  mesh  grid 
for  comparison  with  the  results  obtained  using  a transmitting 
boundary.  The  second  differs  from  the  "standard"  only  by 
using  a transmitting  bottom  boundary.  Lastly,  the  third 
employs  a transmitting  bottom  boundary  and  a reduced  100  by 
10  mesh  grid.  However,  the  mesh  spacing  in  the  third  run  is 
the  same  as  in  the  standard.  The  purpose  of  this  third  run  is 
to  test  whether  the  transmitting  boundary  as  constructed  would 
permit  the  execution  of  a problem  with  a reduced  number  of  mesh 
points  while  faithfully  reproducing  the  results  obtained  with 
a larger  number  of  grid  points. 

The  first  and  second  runs  are  identical  (Figures  13  and 
14)  until  just  after  t = .4l6l  when  the  close-in  waveform  and 
magnitudes  change  slightly.  These  changes  reasonably  coincide 
with  the  earliest  possible  arrival  at  the  surface  (t  = .41 66) 
of  a bottom  reflected  signal  in  Run  #1.  As  would  be  expected 
the  surface  displacements  after  t = .4266  are  more  negative 
in  Run  #2  in  which  reflection  has  been  prevented.  Table  II  shows 
the  surface  vertical  displacements  for  mesh  point  (i  = 1 3 , 
j = 2)  from  Run  #1  (standard,  rigid  bottom  boundary)  in  compari- 
son with  Run  #2  (transmitting  boundary)  in  which  reflection  has 
been  prevented. 
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But,  unlike  the  encouraging  results  obtained  with  Run  #2, 
Run  #3  deviates  from  the  "standard"  at  an  earlier  time  than 
is  expected  from  computing  the  earliest  possible  arrival  time 
of  a bottom-reflected  signal,  with  the  surface  displacements 
becoming  significantly  larger.  This  result  casts  a shadow  of 
uncertainty  over  the  validity  of  the  transmitting  boundary 
being  utilized. 
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LE(CM)  = 0 v ( CM)  = .100E  + D'+ 


Figure  13 


Figure  13.  cont 


Figure  14.  Computations  with  Model  to  Test  Operation  of  the  Transmittin, 
Boundary  - 100x25  Mesh  with  Bottom  Transmitting  Boundary" 


Figure  14.  cont 
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Figure  15.  cent 


Table  II. 


Comparison  of  Surface  Seismic  Motions 
Versus  Transmitting  Bottom  Boundaries 


Rigid 


Surface  Vertical  Displacements 
at  Mesh  Point  (i=13.  j=l) 


0.311 

0.350 

0.375 

0.400 

0.425 

0.450 

0.475 

0.500 


Standard  w/rigid 
Boundary 
(cm) 

100x25  Mesh  with 
Transmit.  Boundary 
(cm) 

100x10  Mesh  with 
Transmit . Boundary 
(cm) 

0.0 

0.0 

0.0 

0.021 

0.021 

0.022 

-0.017 

-0.017 

-0.015 

-0.023 

-0.023 

-0.040 

-0.023 

-0.023 

-0.050 

-0.018 

-0.025 

-0.060 

-0.009 

-0.030 

-0.061 

-0.003 

-0.004 

-0.012 

Simulation  of  a 28  Kiloton  Atmospheric  Burst 

With  the  preceding  analysis  completed,  the  model  is 
next  used  to  simulate  the  airblast-induced  ground  motions  from 
an  atmospheric  burst.  A yield  of  28  kilotons  is  chosen,  with 
detonation  occurring  at  a height  of  burst  of  165  feet.  This 
height  of  burst  is  intermediate  to  those  for  events  in  Table  I 
detonated  at  Frenchman  Flats  and  displaying  an  airshock  precursor. 
Appendix  F gives  the  subroutine  which  calculates  the  airblast 
parameters  on  the  ground  for  this  event. 

The  seismic  reactions  of  four  differenx  ground  media 
are  compiled  in  the  tables  which  follow.  This  compilation 
results  from  the  use  of  the  computer  model  with  rigid  right- 
hand  and  bottom  boundaries.  From  run  to  run,  only  the  ground 
medium  over  which  the  burst  is  assumed  to  occur  is  changed. 

These  four  ground  media  are  described  in  Figure  16. 

The  resulting  maximum  upward  surface  velocities  are 
given  in  Tables  III-V.  These  velocities  were  obtained  by 
differentiation  of  the  vertical  surface  displacements  computed 
with  the  model.  No  maximum  velocities  are  given  for  Run  #4 
which  simulated  the  Frenchman  Flats  test  site  because  all  up- 
ward velocities  were  less  than  0.1  cm/sec  arid  considered 
insignificant.  Failure  of  this  layered  medium  to  produce 
the  enhanced  velocities  indicated  as  possible  in  Appendix  B 
is  apparently  due  to  the  presence  of  the  very  slow  upper  layer 
where  the  seismic  velocity  is  only  slightly  greater  than  the 
sound  speed  in  ambient  air. 
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Figure  16 . Ground  Media  for  Model  Simulation 

of  28  Kiloton  Event 
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Table  III 


Table  IV.  Surface  Velocity  Data  in  Varied  Ground  Media  - Run  #2 


Table  IV. 


Table 


The  stiff,  uniform  ground  medium  used  in  Run  #2 
produced  the  greatest  vertical  velocities.  While  initially- 
appearing  to  conflict  with  Equation  ( 1 ) , this  result  is 
apparently  due  to  the  ability  of  the  stiffer  medium  to  trans- 
mit signals  ahead  of  the  airblast  which  were  created 
"upstream"  at  higher  overpressures,  while  these  same  signals 
are  not  able  to  outrun  the  airblast  in  the  slower,  more 
compressible  medium  of  Run  #1.  Another  factor  is  the 
ability  of  a stiffer  medium  to  more  readily  transmit  the  higher 
frequency  disturbances.  These  higher  frequency  components  may, 
in  part,  explain  the  larger  vertical  velocities. 

Because  the  heating  from  the  fireball  which  any  rising 
dust  receives  is  of  interest  in  creating  the  thermal  layer 
which  determines  the  precursor,  Tables  III-V  also  give 
data  pertinent  to  the  fireball  radiation  emission.  Significant 
amount  of  fireball  heating  can  be  received  by  rising  dust 
only  when  the  difference  between  Tl , the  time  of  occurrence 
of  maximum  upward  surface  velocity,  and  T2 , the  time  of 
airblast  arrival,  is  sizeable.  This  will  only  occur  when  the 
seismic  signal  is  able  to  substantially  outrun  the  airblast. 
Such  substantial  outrunning  is  not  reflected  in  the  Tables 
III-V  because  the  upward  velocities  continue  to  diminish 
with  increasing  ground  range. 

Dust  layers  created  ballistically  would  reach  a maximum 
height  of  only  0.09 6 cm  for  Run  #1,  0.64  cm  for  Run  #2, 

0.16  cm  for  Run  #3.  and  a negligible  height  for  Run  #4. 
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IV.  Conclusions  and  Recommendations 

Conclusions 

Accuracy  requirements,  interference  from  artificially 
reflected  signals,  and  computer  storage  limitations  together 
prevented  the  comparison  of  results  obtained  with  this  model 
(rigid  righthard  and  bottom  mesh  boundaries)  with  the  available 
published  data.  The  data  on-hand  consisted  of  independent, 
elastic  calculations  of  late-time,  airblast-induced  ground 
motions  and  of  late-time  seismic  motion  data  from  actual  test 
events. 

Prevention  of  the  artificially  reflected  signals  and 
reduction  of  the  computer  storage  requirements  would  have 
permitted  this  comparison.  Application  of  the  transmitting 
boundary  condition  constructed  to  accomplish  the  above  two  goals 
yielded  indefinite  results  at  best.  This  uncertainty  in 
whether  the  transmitting  boundary  condition  was  properly  con- 
structed prohibited  its  use  in  validating  the  model  against  the 
independent  data. 

Thus,  the  model  remains  unvalidated,  and  data  derived  from 
it  must  be  viewed  accordingly.  Nonetheless,  it  appears  that 
several  general  conclusions  can  be  made. 

First,  the  surface  velocities  resulting  from  simulation  of 
a 28  kiloton  atmospheric  burst  at  500  feet  height  of  burst  yielded 
a dust  layer  ballistically  reaching  only  0.64  cm  at  its  highest 
point  for  the  stiff  one-layer  ground  medium,  0.096  cm  for  the 
softer  one-layer  medium,  and  a negligible  height  for  the  more 
realistic  four-layer  Frenchman  Flats  medium.  These  heights  are 


\ 
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significantly less  than  the  two  to  three  meter  high  dust  layers 
known  to  exist  prior  to  airshock  arrival.  Therefore,  it  must 
be  concluded  that  as  postulated  (ballistic  rise  only)  the 
airblast  is  not  likely  to  significantly  contribute  to  the 
precursor. 

However,  the  magnitude  of  upward  vertical  velocities 
computed  indicate  that  the  seismic  motions  of  the  ground  surface 
can  be  significant  and  should  be  considered  in  any  modeling  of 
thermal  layer  precursor  generation.  In  particular,  the  seismic 
ground  motions  can  alter  the  velocity  with  which  dust  is  injected 
into  the  air  by  another  causal  mechanism  such  as  the  thermally- 


induced  soil  blow-off  mechanism  researched  by  AFWL.  Only  when 
another  causal  mechanism  can  be  shown  to  occur  sufficiently 
in  advance  of  the  arrival  of  the  seismic  signals  to  preclude  the 
interaction  of  the  two  mechanisms  does  it  appear  that  the  airblast 
effect  can  be  ignored. 

In  addition,  AFWL  modeling  of  the  thermally- induced  soil 
blow-off  indicates  that  the  injection  velocity  of  the  particles 
can  range  from  50  to  200  cm/sec  without  significantly  altering 
the  height  of  the  dust  layer  computed  by  their  model.  At  the 
maximum  injection  velocity  of  200  cm/sec,  the  dust  would 
ballistically  rise  only  20  cm — not  the  2-3  meters  found  expert- 
mentally.  This  indicates  that  thermal  radiation  transport  and 
hydrodynamics  is  the  more  dominant  force  lifting  the  dust  particles. 
Redefining  the  airblast/seismic  model  to  include  these  two  effects 
may  well  result  in  airblast-induced  thermal  layers  more  nearly  in 
agreement  with  experimental  data.  Or  stated  another  way,  the 


key  factor  may  be  to  begin  the  soil  particles  in  an  upward 
motion  whereby  the  thermal  radiation  transport,  hydrodynamics, 
and  Taylcr  instabilities  dominate  to  lift  the  soil  to  the 
significant  heights  of  several  meters.  The  role  of  airblast- 
induced  vertical  motions  in  the  precursor  formation  certainly 
warrants  further  research.  (Prentice,  J. , 197b:  13) 

Finally,  velocities  imparted  to  the  surface  dust  before 
and  after  arrival  of  the  airshock  may  have  considerable  effect 
upon  the  amount  of  dust  swept  up  as  the  airshock  passes.  This 
airborne  dust  would  be  expected  to  have  significant  impact 
upon  the  dynamic  pressure  and  the  erosive  ability  of  the  airshock. 


Recommendations 

Recommendations  for  further  work  with  the  computer  model 
developed  as  an  adjunct  of  this  thesis  research  include: 

1.  the  further  search  of  published  literature 
in  an  effort  to  find  early-time  seismic 
motion  data  with  which  tne  rigid  boundary 
model  may  be  validated; 

2.  the  further  search  of  published  literature 
to  find  one  or  more  suitable,  generalized 
algorithms  for  peak  overpressure  and 

positive  pressure  phase  duration;  such  algorithms 
would  eliminate  the  requirement  of  generating 
a new  surface  pressure  function  for  each  new 
problem  computed; 

3.  the  further  study,  analytical  development, 
and  construction  of  a transmitting  boundary; 
if  successful,  this  work  can  have  wide  appli- 
cations in  this  simulation  as  well  as  other 
models  using  a fixed  finite  difference  mesh; 

4.  refinement  of  the  model  to  confidently  give 
arrival  time  data  of  the  airblast-induced 
seismic  motion. 
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Recommendations  for  further  work  on  the  airblast-induced 


precursor  includes 

1.  the  study  of  experimental  data  to  give 
time  of  occurrence  of  the  thermal  layer 
for  comparison  with  seismic  signal  arrival 
times  in  order  to  determine  whether  the 
airblast-induced  seismic  signals  arrive 
too  late  to  significantly  contribute  to  the 
thermal  layer;  and 

2.  if  the  above  investigation  reveals  that 
the  arrival  of  the  seismic  signals  is  not 
too  late,  incorporation  of  thermal 
radiation  transport  and  hydrodynamics  into 
the  model. 
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Appendix  A 


One-Dimensional  Estimate  of  the  Maximum 


Vertical  Velocity 


Suppose  that  Sx>  the  displacement  of  a particle  in  the 
ground,  can  he  expressed  as  a wave-like  disturbance  propa- 
gating in  the  x-direction  with  velocity  c. 


Sx  = f ( t-c/x) 

The  displacement  S then  satisfies  the  wave  equation. 

A 

This,  of  course,  must  be  the  case  as  stress  within  an  elastic 
solid  is  known  to  be  a generalized  form  of  Hooke’s  Law. 
(Kolsky , 1953*  8) 

This  one-dimensional  stress  can  be  expressed  as  the 
following  (Newmark,  N. , 1962:  C-9) 


stress  q- 

^ X 


(A-l) 


(A-2) 


where  ^ is  the  density  of  the  solid  medium. 
Differentiating  the  function  f yields 


CTx 


a f(t-x/c) 

6 (t-x/c) 


c 

a (t-x/c) 

The  boundary  condition  which  applies  at  the  surface 
x=0  is  that  the  stress  is  equal  to  the  normal  loading  of  the 
pressure  on  the  surface,  or 


t 


' t-x/c, 


b (t-x/c)  x=0 


= -P(t) 


However,  the  particle  velocity  is  simply  the  derivation 
of  the  particle  displacement  with  respect  to  time. 


df  (t-x/c; 
d t 


5 f( t-x/c) 
5 (t-x/c) 


d f ( t-x/c) 
d (t-x/c) 


(A-4) 


X f . 

Substituting  for  — in  Equation  (A-3)  gives 

d (t-x/c) 


c f Vx 


- -P(t) 


■ f 


The  particle  velocity  is  a maximum  when  the  loading 
pressure  is  a maximum.  Finally , 
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Appendix  B 


i 


Effect  of  Ground  Medium  Layering 

on  Stress  (Newmark,  1962s  C-13— C-15) 

Given  the  following  two-layered  ground  mediums 


x=0 ^ 

Layer  1 -p1 , 

Layer  2 f2,  c 2 

x 

f 

where  .p  ^ is  the  density  of  the  ith  layer,  and 

c ^ is  the  seismic  velocity  of  the  ith  layer. 


If  an  elastic  medium  is  assumed  then  the  displacement  of 
a particle  within  this  solid  must  obey  the  wave  equation, 

.2  „ 


a2  s. 


6t‘ 


= -f 


X 


&x‘ 


Solutions  to  the  wave  equation  take  the  form 


Sx  = f (t-x/c)  + g(t+x/c) 

where  f (t-x/c)  represents  a wave  traveling  in  the  positive 
x-direction,  and 

g(t+x/c)  represents  a wavs  traveling  in  the  negative 
x-direction.  - «« 


Now  consider  a disturbance  incident  upon  the  interface 
between  the  two  layers,  namely, 


Sx  = f( t-x/c) 


- - > 
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The  stress  incident  upon  the  interface  is  given  by 


<rx.  = fi°i2  iis 

incident 


= - ~fl°l 


b f 

b (t-x/c' 


(B-l ) 


Let  the  incident  wave  or  disturbance  be  expressed  as  the 
sum  of  reflected  and  transmitted  components. 


= S, 


+ s. 


^incident 


reflected 


transmitted 


or,  alternately, 


f( t-x/c)  = G(t+x/c)  + F( t-x/c) 


From  Equation  (B-l)  the  interface  condition  that  the  sum 
of  reflected  and  transmitted  stresses  must  equal  the  incident 
stress  can  be  applied. 


crx  = <Tx  + cTy 

incident  reflected 


transmitted 


-P  1C1  d f(t-x/cl 
b x 


~P  1c1  c)  G ( t+x/ c ) + -p  2c2 


(B-3) 


R = reflection  coefficient  = constant  in  time,  and 
T = transmission  coefficient  = constant  in  time. 


such  that  Gft+x/c^  = R f(t-x/c1)  and 


F(t-x/c2)  = T fCt-x/cj^) 
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From  (B-2) , . f = Tf  + Rf 


1 - R = T. 


From  ( B— 1 ) and  (B-3),  -m— ■ 

' - -f2c2  -A_P 


:t-x/c. 


(t-x/c2) 


fl°l  -S— 2 

C>  (t-x/c1) 


However, 


c)  (t-x/c 


o ( t-x/c  ^ t 


d (t-x/c^ 


Applying  (B-6)  to  Equation  ( B-5) , 


- f 1C1  4—  = - 2C2  T 


Dividing  by  - ^ ^J, 


+ f ^ R 


1 = 


~f  2C2 

f 1C1 


1 + R = 


? 2°2 


Combining  (B-4)  and  (B--7)  gives 


RIRPPiMP^W 


- ■ + 1 
-f  1C1 


* 


since  G = Rf,  — - = R -~~ 

b t ot 

but 

5 G (t+x/c ) ^ (t+x/c)  £>  f(t-x/c)  ft  (t-x/c) 

— r 

c)  (t+x/c)  d t c>  (t-x/c)  dt 


giving 

-f  1°1  = fl°l  r ? f , - 

d (t+x/c)  b (t-x/c) 

or,  equivilently , 

CTr  = R 0'i  . (B-10) 


Similarly, 

F = Tf 
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Appendix  C 

Derivation  of  Algorithm  for  Decay  of 
Overpressure  with  Time 


To  accurately  apply  the  overpressure  on  the  ground 
surface  requires  that  its  time-dependent  behavior  be 
incorporated.  Accurate  behavior  with  time  is  necessary  to 
couple  the  correct  frequency  components  of  the  airblast  into 
the  ground.  To  this  end,  it  is  undertaken  to  transform  the 
data  contained  in  Glasstone's,  The  Effects  of  Nuclear  Weapons , 
and  presented  in  Figure  1?  into  an  algorithm  which  could  then 
be  incorporated  into  the  computer  model. 

It  is  first  observed  that  the  data  curves  are  very  nearly 
symmetric  about  the  line  indicated  in  the  figure.  Next,  a 
"super  ellipse"  of  the  following  form  and  centered  at  the 


point  (1,1)  is  argued  to  reasonably  fit  the  curves  (Nickel,  1978): 

(l-x)m  + (l-y)m  = 1 


At  the  symmetry  axis  chosen,  x=y,.  Or,  in  this 
application 

P ( t ) /PEAKP  = t/tj 

where  P(t)  = pressure  at  time  t, 

PEAKP  = peak  overpressure, 

t = time  measured  from  arrival  of  shock 

front , and 

t*  = direction  of  positive  overpressure 

p phase. 


(C-l ) 


At  the  symmetry  axis  Equation  (C-l)  becomes 


(l-x)m  = £ 

from  which  results 

m = In  2 

In  (1-x) 

From  Figure  17  data  and  Equation  (C-2),  the  following 
table  results t 


PEAKP 

x = t/t+ 

1-x 

m 

1 

20 

• 35 

.65 

1.60 

100 

.27 

• 73 

2.20 

200 

rH 

CM 

.79 

2.94 

1000 

.12 

00 

CO 

• 

5.42 

After  several  trials,  the  above  table  data  relating 
PEAK P and  m is  chosen  to  be  fit  by  a cubic  equation  in 
(In  PEAKP ) , 

m = 1 + A * In  PEAKP  + B * (In  PEAKP)2 

+ C * (In  PEAKP)3  (C 

Three  equations  are  set  up  using  Equation  (C-3)  and  the 
data  table,  then  solved,  giving 
A = 0.382, 

B - -O.136,  and 

C = 0.025 
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Knowing  m,  Equation  (C-l)  is  solved  for  y yielding 

f .mil/"1 


y = i - 


1 - (l-x)‘ 


or  equivalently, 

P(t) 

PEAKP 


= 1 - 


1 - (1-t/tJ)1 


The  results  of  using  this  curve-fitting  algorithm  are 
given  in  Figure  18  for  the  surface  overpressure  from  a 28 
kiloton  event  at  500  feet  height  of  burst. 


Appendix  D 


Derivation  of  a Transmitting  Boundary 

At  a boundary  it  may  be  desirable  to  pass  or  transport 
seismic  signals  through  the  edge  of  the  mesh  as  though  a 
semi-infinite  expanse  of  ground  material  existed  beyond  that 
edge.  Although  several  conditions  may  exist  upon  which  such 
a "transmitting"  boundary  can  be  constructed,  an  argument 
based  upon  momentum  transfer  is  chosen. 

Consider  the  finite  difference  mesh  near  such  an  edge, 
in  this  case,  the  bottom. 

r 

• =«- 

JMAX-2 
JMAX-1 

JMAX 

i 1+i  1+1 

Assume  that  the  horizontal  level  JMAX-1  is  the  boundary 
across  which  seismic  signals  are  desired  to  be  transmitted 
without  reflection.  The  points  marked  with  the  symbol  "X" 
are  known.  Thus,  the  stresses  rr  and  rr  can  be 
calculated  at  point  P,  centered  at  (i+-J,  JMAX-3/2). 

The  stress  Q~  represents  the  'flux  of  z-momentum  in 
z z 

the  z-direction,  and  rr  the  flux  of  r-momentum  in  the 
z-direction.  A properly  constructed  transmitting  boundary 
would  permit  the  outflow  of  positive  momentum  and  the  inflow 
of  negative  momentum,  while  preventing  the  inflow  of  positive 
momentum  and  outflow  of  negative  momentum  which  result  in  a 
net  increase  in  momentum  flowing  into  the  finite  difference 
mesh. 
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The  z-momentum  density  is  given  by  the  product  -fv  , 
and,  similarly,  the  r-momentum  density  by  , where  V and 

Vr  are  the  velocities  of  the  ground  particles  in  the  z-direction 
and  r-direction,  respectively.  Assuming  the  density  -f  to  be 
constant  near  the  mesh  boundary  of  interest,  the  positiveness 
or  negativeness  of  z-  and  r-momentum  can  be  determined  by 
testing  the  sign  of  the  V and  V velocity  respectively. 

Z 


Then,  consider  the  following* 


dS 


dt 


dz 


dz 


dt 


where  S is  the  displacement  of  the  bround  in  the 
z z-direction. 

Here  ^ is  assumed  to  be  greater  than  or  equal  to  zero. 
The  following  logic  tables  results 


»Sz 

^zz 

P 

indicate 
reflection  ? 

action 

» 0 

=»  0 

Yes 

O’zz 

= 0 

R 

» 0 

«=  0 

No 

&~zz 

R ^~zz 

P 

«=  0 

^ 0 

No 

6~zz 

R =CrzZ 

P 

«=  0 

< 0 

Yes 

^zz 

VO 

II 

o 
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Crz 

p 

indicate 

reflection  ? 

action 

p 

> 0 

Yes 

M 

w 

n 

o 

.=  o 

| 

No 

Orz 

R = °"rz 

P 

=»  0 

No 

C'rz 

R =Cr- 

P 

«=  0 

Yes 

<Jrz 

1 

= o ! 

R 1 

I 


C 0 


Yes 


R 


0 


g:  an  outpjt)  ?n?  76 

r>SIO%  V (102,??,?) ,7<1G0,?$,T>  RNO  7? 

FUSION  ,1TP-M  00)  ,C7P(in0)  ,DZ=(100)  100)  RNP  76 

F>‘?3  0>.'  U(t0n,’l)  jHM  7 9 


Appendix  E 


Listing  of  Computer  Code 


O 


(JL 


✓ : r 

4-«  fi  f-< 

c.  c ». 


jO 

a 

jo  or 

u 

X? 

ai 

t?0 

U.) 

<r. 

cO 

aO 

on 

aD 

an 

a. 

an 

on 

on 

X) 

it  i 

aj 

klD 

F- 

K 

K.  K 

F* 

Fv 

Fv 

!>»■ 

•v 

h. 

K 

»s. 

rs. 

Fs. 

Fv. 

K 

r- 

K. 

K 

K. 

K. 

TV. 

G 

n 

G C 

a 

a 

g 

0 

0. 

G 

G 

X 

D 

G 

a 

0 

a 

G 

a 

t 

0 

0 

c 

G 

ZT 

X 

!!  X 

TT 

A. 

z 

T 

77 

— » 

A 

X 

7 

x 

— » 

X 

X 

—r 

x 

A. 

x 

X 

A 

X 

(Y 

rv 

iV  rv 

rv 

fV 

rv 

rV 

rv 

rv 

rv 

to 

rv 

rv 

r>y 

iV 

rr 

rv 

rv 

Ct 

rv 

rv 

rv 

rv 

Z 

o 

c 

H 

o 

X 

t*  1 

in 

ll! 

H 

cr 

cv 

• 

1- 

<T 

to 

U 

7' 

X 

a 

CM 

U' 

CO 

F— 

C5 

C/' 

*» 

<r 

V 

r^. 

-or 

4—4 

4H 

X 

F— 

llJ 

4- 

o 

C 

C~ 

Y 

ll 

<y 

F- 

</^ 

•k 

cr 

u. 

f- 

C 

C/ 

-j 

r 

<1 

• 

• 

00 

cr 

•- 

-J 

to 

-J 

z 

z 

o 

r 

G' 

V 

C 

M 

L» 

4-« 

o 

rv 

a; 

C 

< 

r 

3 

V 

f* 

<T 

o 

to 

X 

0* 

Y 

O 

X 

o 

^■V 

V 

X 

u 

Q" 

F- 

> 

M 

W 

c 

4H 

>- 

l/) 

c 

a 

X 

G 

o 

<1 

F— 

4-4 

z 

u» 

u 

F-4 

U) 

c 

UJ 

•« 

1/5 

F- 

-J 

— i 

G' 

Q. 

UJ 

O 

o 

ft 

ixT 

<1 

UJ 

-J 

lv 

& 

X 

cr 

17 

-J 

a 

<3 

«J 

4— 

O 

Cl 

o 

tv. 

L 

a 

3 

.V 

> 

F-4 

G 

F— 

W 

X 

o 

-J 

O 

o 

w 

o 

c 

rv' 

o 

u 

3 

U‘ 

<T 

f 

G' 

> 

<T 

X) 

O 

CX 

CVJ 

o 

X 

> 

U 

w 

«s 

3 

c 

•* 

UJ 

> 

n 

F* 

x 

cr 

4- 

c: 

rv 

— J 

r* 

(? 

rv 

u 

7 

♦. 

4— 

O 

X 

T 

CD 

o 

ro 

4— 

CD 

X 

<T 

7 

r> 

4- 

CO 

r 5 

H 

O 

r 

o 

«-r 

j 

»-4 

J— 4 

w 

<3 

X 

u> 

u 

Al 

u. 

X 

— i 

L>. 

G 

U- 

_L 

o 

X 

a 

o 

<r. 

>— 1 

e 

*2 

o 

h 

O 

I 

X 

X 

•^c 

c 

-J 

*k 

F~ 

tv 

>- 

15 

X 

«a 

u* 

• 

cr 

t) 

1- 

X 

4-  - 

m 

F4 

n 

X 

o 

m 

cr 

0- 

3" 

Ul 

rv-' 

td 

r> 

r > 

4- 

3 

4-4 

—k 

ri 

tr> 

< n 

o 

r •» 

4-4 

•* 

n 

a 

rv 

4— 

I«> 

it  i 

•— * 

-r 

rv 

r. 

i<  i 

rv 

4— 

r/, 

r 

F 

U 

r 

17 

C 

Ct' 

C 

x 

r* 

•» 

<3 

00 

3 

u 

«r 

77 

U' 

»> 

C 

o 

• 

• 

x» 

o 

or 

CJ 

o 

i-r 

FH 

1^4 

o 

cr 

o 

•k 

cr 

cr» 

it 

U-' 

7T 

0 

-J 

a 

T 

ro 

rH 

a. 

4— 

•— 

<r 

UJ 

>- 

«* 

IV 

u 

o 

• 

G* 

»-* 

Uj 

<1 

X 

F 

C ) 

o 

r ^ 

F-> 

o 

U i 

e 

X 

e> 

CJ 

a 

G‘ 

•k 

•k 

■>_ 

4— 

X 

M 

1—4 

<z 

X. 

UJ 

Z5 

o 

m 

UJ 

U 

UJ 

O 

-J 

Ul 

o 

4-4 

K 

37 

<r 

h-4 

<x 

•» 

c 

C' 

c 

U 

C 

u 

c 

O 

Ul 

o 

(/J 

t/j 

•k 

u 

07 

u 

ir\ 

( 

O 

r/c 

f 

f ' 

A. 

r 

T 

♦ — 4 

y 

fv 

O 

L 1 

-J 

o 

r ' 

r> 

4-* 

»- 

f/i 

4-1 

O 

rv 

ll 

1-4 

u 

_i 

L 

1 

c 

t - 

F- 

c. 

O 

UJ 

X 

1 

li 

V 

4— 

h 

r 3 

4 — 

1 

n 

X 

o 

< 

4- 

«_ 

L. 

l< 

I- 

f 1 

c. . 

F 

•k 

F-* 

"T 

w 

n 

r*. 

▼J 

a 

«r r 

1-4 

X 

•FT 

F> 

— J 

•s 

cr 

rv 

x 

^ r 

c ■ 

rr  ^ 

o 

■/5 

r-J 

_l 

CJ 

!*• 

7} 

<1 

r 

T 

TT* 

U' 

r 

o 

O 

o 

cd 

r ^ 

FH 

L 

It 

—1 

t 

4-— 

o 

4— 

•> 

3 

— i 

UJ 

n 

7/ 

rv 

x-« 

X 4 - 

c 

4— 

» 

V 

c 

i 

3. 

tl 

a 

u 

F 4 

i. 

li 

• 

i .• 

w 

_J 

c 

XT’ 

111  ’ 

- 

t ) 

i: 

Cl 

c, 

a, 

U 

L 

o 

F4 

> 

K 

rv 

k.j 

NJ 

V 

X _J 

V * 

V.’ 

:i 

' -1 

4-1 

I 

(J 

k 

r~ 

»- 

> 

i . 

_ 

h 

•* 

v n 

O 

X 

_J 

»• 

i— r 

a 

»k 

CD 

•k 

►« 

a 

«7 

X 

^ K 

.w* 

c. 

l 

( 

4 1 

*r-4 

T- 

C' 

; 

r * 

wj 

m 

f 

C) 

a 

X 

r- 

i"- 

4—4 

F 1 

T' 

X 

s 

C 

w 

X 

s 

T» 

♦— 1 

Ft 

F- 

>- 

Ct 

l 

O' 

C 

ll 

C 

CO 

to 

— j uj 

n 

4-4 

X 

Ik 

4- 

; 

H 

h* 

-V 

*. 

A. 

g 

j. 

#k 

h 

•a. 

1 

c* 

r 

V_ 

F 

u 

ll 

-J  u 

c 

• ./ 

4-< 

C 

3 

*-k 

x 

1 

c 

1\ 

V 

1;: 

UJ 

r- 

/. 

l F 

'-l 

i 

> 

4-  4 

i/j 

»•  • 

. _r 

4 -4 

w- 

0 

• H 

F-* 

fH 

1: : I* 

l i 

r 

“0 

0' 

X 

r • 

CL 

c 

fr 

o 

w 

V. 

c. 

C 

<i  r— i 

a 

u 

i. 

7,. 

U 

u 

a 

Ll 

Li 

o o o 
ecu 
c c < 

c c e, 

f f f f t G t f f ' C ' f f t'  C‘  f 

c c C f C ft  c c c c t c c t 


«•  » t>  ♦ V 1 


4*"1 


x-  r r 

r r.  n 

G-  r r> 


17-R0=  3P0MMH  OISTANC-  rR04  G'OJND  7rR0  AT  WHICH  SHOCK 


a_>  cn  to  cn 

K-  4s^  X K* 


aooQ.  qoqc  r.  cr*  a 

777777^77777 

rv  rv  rv  rv  cv  iy  rv  rv  rv  rv  rv  rv' 


cn  on  co  cn  m co  a a : cn  nr 

K.  K.  K K.  N.  N.  K K K.  K 

aoancaoQau 

77?7r27?77 

ivrvivft'ivo'iyo'o  (V 


Ji  <fi  m tf/  «j  <r,  hj  <o  v£>  <d  co  a d 

K K K.  h ^ ^ K K_  N. 

or.  Goari'-onr.  acioc. 

. ’■£  «-  . Li*  * " ” ^ * — /'  ^ 

o rv  iv  cy  (V  cv  <t>  <v  rv  d iv  rv  rv  rr' 


3 

3 

3 

-J 

X 

X 

X 

X 

II 

fO 

d 

2 

• 

7> 

ro 

•v 

O 

CVJ 

I- 

rv 

rv 

rv 

r. 

3' 

4-1 

IxJ 

O 

• 

• 

• 

• 

* 

• 

o 

rv 

• 

(/) 

L! 

<t 

#v 

r. 

rv 

rv 

W 

r 

M 

V 

ro 

K 

V 

▼h 

X 

• 

• 

• 

r 

II 

X 

Lx 

(S 

•» 

3 

• 

• 

CV 

<>~V 

CV 

CM 

(V- 

r 

rv 

•v 

2 

cr 

O 

r4 

K 

K 

tH 

ro 

r-w 

X 

X 

• 

2 

r 

Li 

• 

L 

r 

Li 

r 

a 

r 

<1 

o 

a 

• 

■ 

•s 

CO 

rv  O 

rv  CD 

rv  O 

2 

3* 

lit 

«X 

rv 

X 

cn 

o 

■ 

• 

2 

rl 

2 

rl 

2 

ri 

rv  3 

CO 

> 

‘t 

o 

\Y  CY 

a 

a 

a 

LJ 

X 

4H 

• 

Ll 

LD 

rv 

r- 

r» 

•v 

• 

J 

C* 

CM 

V 

^“v 

• 

• 

• 

• 

CD 

r 

• 

rv 

Q: 

r-H 

LJ 

<c 

IV 

(No 

CV 

CV 

CV 

rH 

r» 

ro 

d 

LJ 

(Y 

«J 

L 

Hi- 

It 

4- 

1* 

K 

II 

li 

b< 

ro 

• 

0v 

<3 

O 

ll! 

c 

V 

# 

r* 

4—i 

v.3 

V 

m 

L. 

t: 

r: 

rn 

X 

ro 

7 

ro 

x 

X 

• 

rv 

L 

a> 

iY 

01 

O 

3 

a- 

< 5 

V 

a 

•i 

c. 

r 

* rv 

.^W 

o 

HH 

(Y 

V 

\ 

4- 

V 

V 

V 

•4 

II 

2 

• 

X 

c: 

L 

C/t 

C 

3 

co 

3 

u 

T 

cn 

3 

II 

II 

A 

tn 

> 

X 

LI 

ll> 

r •» 

It 1 

f o 

ll  1 

r *> 

it 

t > 

4— 

x 

cr 

1— 

X— 

V. 

*r 

V 

X 

x 

x 

V 

3 

o 

Lx. 

l_^ 

_i 

a 

«a 

> 

a. 

> 

X 

> 

> 

C. 

<J 

a 

W 

Of 

7 

h-' 

O 

C7 

LT 

C 

O 

a 

Lr 

03 

C 

X 

1 

O 

LI 

a 

3T 

O 

w 

w 

w 

W 

X 

4-1 

4- 

L. 

►-! t 

o 

_J 

>- 

>- 

3- 

> 

Cl 

*— 

■ 

L) 

x 

1— 

4— 

i— 

W> 

l_ 

4-< 

cr 

TT 

tn 

D 

it : 

II 

4-' 

li 

4-H 

n 

»-1 

II 

4-1 

rv 

2 

n' 

C'l 

“T 

"S' 

<v 

f/i 

n 

IA 

f/> 

II 

#/• 

/A 

rv 

n 

• 

io 

d 

n 

rr 

d 

X 

•^r 

X 

«-l 

rv' 

♦ 

r 

n 

l 

r 

l: 

c 

l 

r 

ll 

r 

t 

L 

r 

x 

H 

**R 

X 

O 

L 1 

O' 

CO 

3 

n 

T 

o 

r 

n 

r 

O 

4- 

vD 

4H 

LJ 

CO 

»-♦ 

X 

n 

<1 

Lx’ 

d 

<r 

d 

<3 

Lu 

L 

CO 

0V 

C/J 

H 

O 

•J 

-j 

a 

a 

X 

2. 

r- 

LJ 

z 

Li  1 

c 

a 

cn 

<1 

r* 

r> 

rv 

rv 

rv 

rv 

<T 

«r 

2 

C 

r> 

L 

• > 

7 

a 

r- 

• 

a 

2 

rv 

2 

f r 

• 

5 ' 

if 

II 

a 

7. 

(7 

j) 

K- 

Q 

(Y 

rv 

LI 

rv 

f> 

rv 

L J 

r> 

— i 

o 

O 

r» 

y- 

a 

a 

U j 

re 

>- 

ro 

L 

r 

>■ 

X 

rv 

Q 

ll 

4— 

c 

t 

cn 

cr 

71 

Lit 

i- 

>- 

• 

O 

d 

r 

O 

V 

• 

05 

d 

• 

X 

41 

«■ 

C 

X 

c. 

c 

cv 

n 

c 

7 

« 

r 

7 

CV 

c ; 

-J 

L 

U! 

LI 

L 

4- 

K 

4—i 

X 

-J 

rH 

CY 

L 

-J 

u 

y\ 

2 

«u. 

a 

X 

M 

4-H 

c~ 

n 

tn 

IX 

14' 

li 

n 

-J 

l. 

t t 

li 

r 

a 

L 

4 

X 

ti 

4— » 

4- 

r* 

L 

T 

«• 

r>" 

d 

-J 

r\ 

r 

4-* 

r 

r> 

r 

a 

r- 

r- 

<- 

tn 

4- 

x 

r 

o 

C 

* 

L 

r\ 

l-H 

• 

75 

r. 

• 

7 

l-i 

• 

T' 

c 

2 

1- 

4-1 

V- 

1— 

CT 

Li 

0v 

< > 

> 

O 

3 

fO 

5* 

r> 

X 

</' 

r 

4— 

II 

Ci 

• 

L. 

«. . 

*7 

& . 

3 

it 

n 

1 » 

(/ . 

( 

r 

4— 

li . 

II 

n 

X 

r-> 

c 

cf 

M 

r. 

r 

11 

rv 

II 

r 

X 

II 

a 

X 

4 

L> 

Ui 

C 

o 

fj 

cn 

4 

X 

1-1 

•v 

Cl 

T 

7 

X 

C 

x 

u. 

r*« 

4-  < 

i . 

r° 

a 

a' 

n 

75 

h- 

I- 

X 

li 

c*r 

O 

■? 

c 

r.' 

•v 

< 

i0' 

4-H 

L. 

Li 

* 

it 

i. 

5 

<1 

.1 

a 

CV 

-J 

r— » 

CM 

1 

7^ 

> 

• 

■j 

-• 

r 

rsi 

to 

c 

L 

« . 

T 

_J 

(I 

*- 

L 

L. 

- 

c 

3" 

• 

r 

( 

vi 

( — 

:< 

»- 

£ 

n 

r "* 

C > 

/ 

r 

C ^ 

1- 

d 

rv 

(/* 

•v 

1— 

r- 

li 

r> 

r > 

•V 

i ~ 

•» 

x 

r 

* 

r. 

V 

ri 

T 

T-* 

7 

(V 

x 

r 

U 

O 

(/ 

. 

f 

, 

fO 

r* 

f 1 

Cl 

. , 

c * 

._T 

v 

f 

cu 

4 ■> 

c. 

vl. 

^•i 

u 

LJ 

G 

2 

l: 

X 

2 

c 

: o« 

i 2 

U * 

2 

\ 

Cl' 

2 

\ 

C.) 

2 

\ 

C.! 

2 

v 

o 

2 

r 

2 

✓ 

C 

X 

w * . 

s-' 

lv 

O 

le 

O' 

r-  • 

■v-' 

0 

*r-< 

/> 

r-H 

V— 

a 

r H 

*— 

i— 

>- 

4- 

1- 

h- 

--9 

h- 

L. 

h- 

L» 

4— 

L . 

4" 

L i 

4-- 

r- 

r - 

C i 

l 

<4 

L 

1 

•-* 

i 

•-C  . 

C ^ 

t 

»* 

-.. 

1 

«-v 

L. 

4 - 

1 

•vl 

L 

I 

;■  v)  r > — ~ r n : ir 

Ci.’  «3  M .1.  t-l  i r-H  y >-■ 

o r.  o : -i  ■ o r-  o rv 

L i/J  L Ir  L b Cl  li  Li 


v v r > 7 5 

O file  O •— * cl 

^ 'V  o w > o 

LI  L L 


^ 2 > S X 

_ .3  »-i  It 

>■  o w o o 

L L L li 


L 

4 1 V hv 

V 3 f— 


pro-  TIM-  0-  SHOCK  A PP  T 7 A L U ( ?LA3ST0NE  OR  0NA3?8iF) 


AD-A063  934 


UNCLASSIFIED 


AIR  FORCE  INST  OF  TECH  WRI6HT-PATTERS0N  AFB  OHIO  SCH—EIC  F/6  10/3 
THE  NUCLEAR  AIR-SHOCK  PRECURSOR!  A STUDY  OF  THE  CONTRIBUTION  OF—ETC<U> 
DEC  78  R N PRICE 

AFIT/6NE/PH/78D-21  NL 


R7ERC  IS  330'JVD  R 4 N r E flT  WHIS4  IT  IS  DESIRES  TO  START  D RO  3 L E M ! RNJP  TE 


cn 

as 

cr» 

cO 

c c 

cf 

cc. 

a) 

CO 

on 

or 

CO 

on 

C f) 

TO 

cr> 

co 

cn 

CO 

cO 

•r> 

ai 

an 

CO 

00 

X! 

CO 

CL' 

cn 

CD 

<r> 

■n 

CO 

K. 

K 

h- 

K 

K 

N. 

N. 

k 

tv 

N. 

K. 

K- 

K 

h*. 

►v 

K 

Kw 

K. 

K 

rv. 

rv 

Vv. 

Vv- 

N- 

a 

0 

a 

n 

t 

n 

C 

D 

a 

0 

0 

0 

Cl 

a 

0 

0 

a. 

a 

a 

Q 

Q 

0 

n 

Q. 

C 

a 

a 

fl 

C 

a 

0 

Q 

0 

x 

X 

* 

ZZ 

—* 

x 

2! 

x 

I? 

Z9 

■7* 

X 

x 

-»r 

X 

z 

X 

2 

«*_ 

X 

x 

z 

X 

X 

3T 

AV 

rv 

(Y  !Y 

rv 

A 

fV' 

rv 

rv 

rv 

f*' 

rV 

fr 

CV 

rv  rv 

O' 

rv 

rv 

rv 

rv' 

rv  o' 

rv 

O' 

rv 

rv 

rv 

rv  rv 

rv 

rv 

rv 

o 


X 

x» 

o 

o 

UJ 

a 

O 

c 

h- 

v- 

V 

X 

h- 

o' 

H 

i-C 

cr 

UJ 

00 

>- 

oo 

>- 

x: 

1 — 

h- 

• • 

ar 

00 

c 

o 

cr 

o 

«-r 

o 

c 

X 

o 

u 

►- 

ZZL 

IV 

z 

O' 

M 

►H 

h- 

h~ 

• * 

0^ 

M 

o 

^V 

a 

o 

a 

rv  n 

ft 

1- 

h 

<1 

<4 

X 

X 

u 

or 

X 

h* 

Cl 

a 

< 

CJ 

CJ 

Ui 

e. 

T 

0 

> 

o 

X 

CO 

co 

U; 

U' 

X 

X 

u 

ft— ( 

h- 

rv/ 

CO 

O 

UI 

>- 

UJ 

r 

0 ' 

Qr 

c 

o 

r 

a' 

UJ 

o 

UJ 

H 

UJ 

U i 

ui 

rv 

Or 

M 

VH 

M 

M 

o 

C/) 

00 

>- 

2L 

or 

rv 

u_ 

CJl 

a 

C- 

Ci 

V— 

t— 

UJ 

a 

zr 

<f 

_J 

c 

X 

<1 

O: 

Q. 

X 

V-i 

X 

o 

o 

o 

o 

1 

1 

CD 

o 

o 

X 

u. 

ZK 

ZD 

O' 

UJ 

-J 

u 

Of 

U 

I— 

_J 

00 

c 

o 

o 

.if  ^ 

UJ 

a 

z 

h 

sr 

X 

cl 

UI 

a 

> 

UJ 

>- 

o 

<4 

U: 

L— 

h— 

i- 

a' 

“0 

& 

a 

tu 

X 

2. 

fv 

-J 

o 

> 

•J 

V-i 

»_) 

ro 

V-- 

X 

ZD 

X 

rv 

p 

V-t 

rv 

Cj 

u 

o 

ftH 

ft- 

o 

M 

-1 

«u 

-J 

X 

r-H 

o 

o 

O 

c 

Cl 

IH 

V-C 

Ci 

h-C 

o 

a 

ft— 

a 

o 

x 

r 

l 

l — 

fO 

cr 

v-t 

»- 

H 

u 

u 

*— i 

UJ 

rv 

<r 

00 

— 1 

-J 

-J 

CV 

fH 

X 

7 

X 

ro 

00 

_i 

X 

zr 

u 

f.o 

c 

u; 

u. 

LJ 

Lh 

_l 

t-H 

_J 

u. 

00 

M 

7 

V— 

t— 

<T 

o 

— i 

o 

u 

rr 

r> 

CL 

h- 

O 

r-t 

o 

UJ 

> 

X 

t-H 

A' 

fK 

zr 

zr 

o 

v—t 

ci 

hH 

Ci 

O 

Of 

u 

d 

cr 

O 

i — ■ 

IA 

o : 

O 

{/. 

c 

(0 

•)- 

G 

>- 

a 

u 

it_ 

M 

V- 

»-< 

h 

11 

ft— 

a 

f/> 

rv 

27 

O 

fA 

y 

v-l 

V 

<r 

5 

X 

u. 

r> 

VH 

III 

rv 

tit 

ft— ( 

cr 

fV 

UJ 

0 

<cr 

T 

m 

T 

...i 

<i 

X 

«3 

rv 

fA 

c 

fA 

U. 

U: 

O' 

00 

CJ 

CO 

ft— 

(V 

X 

n 

< i 

»v 

V— 

u_ 

i* 

C_ 

V- 

v. 

1 

IL 

> 

U. 

U_ 

t_t 

u 

La 

u 

C2 

u_ 

o 

M 

*> 

C_4 

ft“1 

N 

2 

a 

X 

CO 

UJ 

o 

X 

Q 

CJ 

-0 

<1 

<4 

>> 

«SJ[ 

C 

> 

a 

c 

ZT 

• 

X 

a 

ft-H 

1- 

2T 

a. 

o 

M 

M 

ZJ 

il 

-J 

II 

H 

HH 

a 

-J 

M 

a 

ft— c 

•s 

_i 

ft-4 

O 

a 

u. 

X 

O 

K- 

o 

U^ 

"0 

V 

CD 

r< 

CO 

Q 

n 

X 

X 

X 

X 

a 

ht 

2 

• 

('• 

i 

~r 

r v 

fA 

1— 

li 

rv 

i— ( 

V' 

r> 

f/t 

(/• 

M 

r/> 

1-4 

(A 

1-4 

ft — 4 

UJ 

h 1 

<rr 

u 

V- 

<r 

Li* 

r- 

U 

n' 

h- 

t\' 

K-C 

M 

if 

ll 

7 

ft— 

— > 

~T 

c •» 

it 

It 

It 

It 

o> 

«rt 

T“ 

l— 

» 

*7 

-v 

ii 

lit 

r. 

r. 

) 

X 

7 

ft— i 

'W 

r'l 

X 

1 

*a 

c\ 

rv 

i— < 

f A 

rt 

rt 

n 

1 

*—> 

n 

CD 

ro 

< »v 

ft— 

ft— * 

rv 

o 

•> 

li 

ft— i 

W 

1 - 

u 

(1 

IT 

c 

X 

U 

L. 

i 

h- 

c 

T_ 

r 

z~ 

U 

•— 

1- 

u 

Z 

ft— 

ft- 

<1 

0) 

cr 

O'* 

(/' 

(/) 

m 

O 

h- 

C^ 

lL* 

»— 

k 

u- 

1 — 

C 

o 

C'’ 

r3 

to 

— - 

o 

L ! 

ft— 

CJ 

<4 

o 

«■; 

c 

UI 

U. 

U 

U_ 

«J 

• • 

z 

Uj 

o 

zz 

aL 

ZZ. 

-J 

f— c 

HC 

VH 

7 

•n 

ft—1 

”T» 

o 

X 

V- 

tL 

)-< 

rH 

v-i 

Cj 

00 

H 

X* 

<r^ 

y-t 

tH 

hH 

hi 

a 

a 

u 

a 

a 

r/ 

h- 

Z 

ft— 

O'. 

a 

u 

ZD 

2 

U! 

C: 

L— 

k— 

v~ 

V- 

c 

h- 

a 

U. 

rv 

Cl 

o 

■^r 

o 

a 

HC 

V 

X 

o 

Ll' 

ft- 

z 

a 

L 

X 

rr 

a 

it 

• • 

G 

a 

_J 

tV 

Q 

o 

L : 

C 

r— 

a 

r/i 

CO 

c*. 

u 

) 

D 

~T 

Q 

VH 

L 

CL 

Uj 

li. 

b_ 

UJ 

of 

<4 

~T 

<v 

U' 

e 

o 

H 

-J 

h- 

-J 

C» 

ft-c 

a 

o 

o 

r> 

r* 

a' 

a 

Cu 

a 

a 

3 

o 

X 

L2! 

K 

00 

l~ 

u] 

h- 

a 

zz 

IT 

H 

X 

V-c 

a 

a 

u 

.J 

ar  o 

00 

X 

c^' 

-J 

o 

o 

a 

o 

1 — 

u. 

l>0 

t-H 

z: 

O 

h 

O 

UJ 

zz 

I'J 

ZT 

o 

c z 

tr 

•* 

u 

0 

c. 

•— 

L- 

L 

C. 

L 

LT 

o 

0 

4ft' 

J 

vL 

l 

c 

C 

i 

c 

>* 

C 

a 

u 

c 

z 

c 

hv 

x 

UJ 

C 

X 

a 

a 

CL 

a 

V 

r 

tu 

h 

H 

X 

z. 

a 

c 

a 

U ! 

u 

a 

CI 

co 

X 

ft-H 

zz. 

a 

*-» 

<v 

X 

Cj 

Ti 

Z9 

rr 

VH 

r 

o 

o 

O 

o 

CT 

r t 

H 

ft-T 

ft— 

ft— 

c 

<3 

_J 

r 

ft-* 

rv' 

V-’ 

u» 

it 

It 

u • 

tj  • 

» 

r\' 

r\’ 

n 

*7 

r*' 

T* 

f.r 

ro 

c" 

( ° 

rr 

rr 

O 

V 

rv 

O' 

-!► 

rl 

u. 

h- 

Of 

IV 

Ct 

_J 

T 

O 

O 

n 

0 

Ll! 

Li 

_l 

I- 

-J 

y— 

c. 

c 

Ci 

Li 

cr 

o 

— > 

ro 

ft- 

rr 

er 

•-1 

cj 

«:T 

_ J 

u 

• — < 

L 

X 

—v 

cr 

f A 

> 

f/0 

V 

y~t 

VH 

n 

M 

r, 

ft— i 

rv 

ft- 

-0 

c . 

•v 

X 

h 

r» 

—j 

ti, 

1 > 

r. 

J 

wr_ 

_j 

i 

l A 

:r 

f A 

r/i 

ti* 

1(1 

U' 

Jr 

X 

r* 

n 

r* 

O 

C > 

<T 

ti  1 

It* 

• • 

X 

c-r 

<4 

•rr 

n 

r. 

HI 

t — 

t-c  1- 

f 0 

»"C 

Li  ' 

CO 

ft-*  Uj 

•• 

► 

n. 

•i 

o 

2 

X 

X 

x 

7 

CL 

Z" 

X 

r— 

to 

li  . 

tiJ 

3t" 

7 

o 

"I) 

r 

lii 

u* 

b 

cr 

r 

co 

— J 

r,-' 

n 

cv 

Ci 

-J 

V-l 

►-I 

; — 

u 

* 

• » 

T; 

•« 

M 

h-i 

o 

o 

h-- 

(V 

ft-4 

ft— 

< "ft 

o 

u. 

r* 

rx 

<4 

<1 

r 

C.-) 

c 

X 

r- 

Cl 

X 

K 

II 

O' 

it 

r. 

£ 

C 

»-I 

ft- 

i: 

l.  i 

o 

“ 

<. 

*- 

•» 

H 

U., 

«4 

tf , 

h- 

c? 

C 

cn 

O 

o 

CJ 

'■J 

vl 

a 

ll 

ft  4 

o 

ii 

H 

h 

7 

1 

x 

— ■ 

— - 

f 

■ 

— 

n 

\-  ? 

X 

lL 

.J 

_J 

> -4 

r.n 

ft  • 

rr 

c 

fl 

ii 

“S 

r. 

T" 

z~ 

z~ 

I* 

iT 

l-H 

O 

o 

It 

n 

li 

n 

II 

II 

H- 

V- 

i- 

1- 

f'' 

ft- 

O 

c 

— D 

CD 

6 

tU 

f] 

CJ 

Ci 

1- 

l< 

i i 

UI 

il 

a 

i J 

o 

go 

L 

► 

•» 

•x 

j - 

h 

O 

CD 

t— 

h- 

X 

*T» 

HH 

► i 

u 

Z. 

*- 

> 

II 

II 

X 

.H 

1 ^ 

<T 

«u 

o 

Cj 

•• 

•* 

0O 

00 

UI 

ii 

<1 

«u 

<L 

cr 

UL 

a 

-J 

-J 

-.J 

^1 

hH 

h-t 

H-4 

ft  < 

rv 

i 

> 

L 

z. 

U- 

2- 

u 

u 

u 

u 

W- 

0 

o 

_J 

<4: 

n-i 

CJ 

n 

O 

r. 

1-4 

h- 

"0 

V-' 

*0 

LL' 

r' 

r 

> 

M 

C' 

-> 

IV 

c*  o c cj  c:  o c e c c c o c a c o <j  o c.  c e cc.cc  c cj  c c.  c:  o c o c c^ 


/ 


ao 

cO 

co 

on 

cc: 

CO 

X* 

T) 

<£ 

cn 

CO 

cO 

r^. 

K. 

r>~ 

K- 

K- 

K 

N«. 

K. 

a 

0 

a 

0 

0 

C 

a 

n 

Q. 

a 

Q 

Ol 

0 

7* 

7 

2T 

TT 

r? 

7 

IT 

2T 

X 

2? 

2T 

7 

ft' 

rv  rv' 

Q' 

rv 

ft 

rv 

tv 

<v 

rv 

(V 

Tr' 

ft' 

3T 

:r 

O 

• 

• 

• 

-J 

o 

-j 

c 

L 

ui 

I*  i 

L* 

cd 

t- 

CP 

1— 

n 

ft' 

r 

cr 

r 

C 

n 

o 

T 

r 

x 

21 

21 

X 

H4 

hH 

■ 

— 1 

M 

; 

X 

X 

1- 

LlJ 

h~ 

LJ 

*H 

r* 

X 

►H 

CO 

X 

jC 

H- 

X 

»- 

x 

a 

V 

a 

00 

d 

LJ 

CO 

d 

L- 

u_ 

2L 

c. 

lx 

*- 

c 

♦-4 

w 

_J 

lx 

X 

-J 

L 

X 

O 

(/■' 

O 

on 

L.J 

d 

LJ 

Lw 

d 

lx 

o 

li- 

-t.. 

<J 

L 

n 

ft 

CJ 

fV 

lx 

LlI 

f/1 

li- 

LJ 

00 

<V 

l— 

o 

(V 

1— 

c 

IxJ 

21 

LJ 

lx 

21 

lx 

f~ 

M 

liJ 

l- 

M 

LJ 

X 

O 

X 

O 

f-H 

tu 

X 

(VJ 

V-H 

LJ 

X 

T 

LJ 

T 

liJ 

O' 

V- 

4 

CY' 

W- 

ii. 

rv 

III 

(V 

V- 

*- 

Li 

^V. 

CJ 

X 

T' 

Itl 

r-r 

C 

X 

r 

< 

C 

> 

-J 

U 

<1 

O 

o 

-J 

(X 

<r 

IL. 

— 1 

x 

* 

rl 

li. 

-J 

X 

o 

<X> 

ci 

O 

N- 

d 

rr 

CT1  C 

3 

J 

h- 

' Ct' 

CP 

c» 

X 

o 

o 

CJ 

2*1 

-J 

1 

i: 

»H 

o 

C_- 

o 

~> 

o 

3 

II 

L» 

hH 

v~l 

-r-< 

V-i 

II 

LU 

»- 

O 

cr 

• 

or 

»- 

h- 

X 

C2 

00 

o 

k 

Li 

a' 

a 

c: 

h- 

»“ 

a' 

u 

r 

0 

Li 

*- 

L 

2. 

l 

* 

-J 

»-» 

x 

*~V 

>r 

-J 

HH 

X 

* 

ll 

Cl £ 

1 

<r 

1 

9 

L 

ry 

1 

».T 

* 

Ij 

n 

I 

-J 

fO 

4- 

If 1 

n 

1 

-J 

1* 

(V 

1 

h 

1 

> 

O' 

1 

▼H 

r 

X 

• ' 

X 

J 

X 

V 

CJ 

V 

• • 

V 

V 

X 

c 

7T 

<r 

c 

<r 

4 

CN 

T“ 

O 

17 

c 

«.*T 

c 

X 

►H 

T* 

.7 

r-i 

2~y 

V" 

V-* 

X~ 

►“ 

X 

r 

*— 

sy 

l- 

h-' 

r— 

Cl 

• 

20 

• 

Ll 

11 

• 

~ 

iS 

• 

22 

• 

t; 

lx 

h- 

L 

l*J 

LJ 

•-j 

Li 

I — 

•^4 

H- 

O 

r'-' 

oo 

CJ 

l9* 

X 

m 

ir 

a 

r/ 

0 

00 

• 

• 

«r. ; 

9 

• 

• 

CL. 

a 

CL 

rr* 

— 1 

( > 

C' 

l J 

c 

c.. 

• 

C 

Cl 

• 

c/ 

11 

II 

(V 

J 

o 

r ■ 

• 

u 

W. 

Xt 

lx 

v-y 

L. 

*-4 

L) 

*- 

h 

U. 

LJ 

H- 

X 

f- 

X 

0 

o 

>- 

1— 

H- 

>- 

X 

’M. 

C> 

«_ 

u 

4. 

c^. 

X 

* «. 

«. 

c_. 

2* 

l-H 

> 

-J 

L. 

L 

J 

►— i 

- » 

2 

-J 

(J 

•jf. 

V_y 

X 

> 

u 

u. 

X 

e 

L. 

a. 

CJ 

L 

c£ 

L 

O 

L 

L. 

C_ 

L. 

L 

HH 

u 

1 — i 

—I 

-J 

r-i 

U. 

1 — i 

u 

li 

t>> 

CL 

*c 

CJ 

r 

C 

i 

c. 

C- 

c 

▼-* 

rl 

alaicf  co®u®B*ir!cti(n©ir 

KKK. 

a apoanoaoD.  Qoca 

Z7Z7727Z7227  7.  Z 
rv  <y  o o-  .v  y r>  ry  iv  (i  ft  tv  tv  tv 


o 

T-t 


o 

cr 

O 

\ 

i 

1 

LJ 

o 

> 

o 

ln 

«-! 

<r 

•* 

«v 

y 

«-< 

V 

4 

X 

L 

r 

d 

d 

-r 

d 

d 

X 

*tT 

c 

3“ 

r/ 

o 

V 

C^ 

/ 

7* 

X 

”7 

li* 

_J 

27 

7 

“0 

ft 

— J 

ci. 

2 

li. 

• 

X 

L 

2. 

ii. 

• 

L 

•n 

c 

<a 

L' 

c 

a 

O 

L 

X 

(N 

f\ 

II 

00 

_ 1 

It 

0^ 

CJ 

II 

II 

It 

II 

• 

li 

n 

It 

• 

1-  1 

,-v 

~v 

...» 

n 

(1 

*> 

Li- 

D 

r*«* 

Lj 

a 

CJ 

CL 

•» 

“3 

"O 

22' 

L. 

CJ 

”0 

'7 

2J 

Ll. 

CJ 

rH 

CJ 

r— 1 

r. 

X 

X 

CJ 

OH 

•s 

ZL 

X 

u 

T-l 

▼ 1 

w 

H 

r-i 

i — > 

**' 

r—« 

r 

i-  ■ 

«-i 

lx 

T-f 

r 

w 

1- 

-J 

ftrH 

T~. 

w 

r— 

-J 

X 

d 

x» 

Cj 

.~~ 

w 

X 

^2 

O 

/2 

d 

o 

O 

-J 

T 

T 

o 

ll 

C> 

u 

T 

X 

o 

L 

-J 

o 

CT 

u' 

o 

r-i 

r 

f- 

o 

c 

£V 

89 


CL. 

cn  u) 

ac. 

cn 

ao 

cn 

~c 

lO 

u> 

cr 

on 

CO 

CD 

aO 

aO 

uo 

CO 

09 

ax 

00 

GO 

*Ci  CO 

cn 

<o 

V 

K-  N- 

K*. 

N- 

r- 

N 

K- 

*. 

K. 

K. 

r*- 

K K 

iO 

Ct  0 

r* 

0 

0 

D 

a 

0 

a 

a 

a 

a 

2 

2 

2 

2 

2 

7* 

2 

2 

2_ 

a x 

2 

7» 

cv 

2 2 

2 

2 

2 

2 

zz 

2 

2 

cx 

2 

"Z 

X 

X 

X 

X 

X 

X 

X 

X 

X 

2 X 

X 

X 

X 

o 

<V  & 

O 

Ci 

(V 

f* 

rv 

rv 

rv 

tv 

tv 

a 

CO 

€0 

cr 

40 

10 

4 0 

CO 

4 r 

10 

<v  to 

to 

40 

lO 

Jt 


cv 

II 


X 

c 

I 

X 


r 

c 

O 


o 

C: 

b. 

V 


O 

h 

X 
o cr 
l*  7“ 


00 

bJ 


< 

!> 

X 
r/^ 
L ' 
T 


O 

CL 

UJ 

u 


2 

Ul 

X 

Ul 

w- 

cl 

in 

to 


in 


UJ 

2 


bJ 

2 


IV 

O 


1- 

i-i 

V 

bj 

Ii 

a 

c 

r 

<r 

C 

X 

T 

o 

ncf  rt 

O 

V 

a' 

<T 

o 

K 

1- 

UJ 

-J 

y*** 

X 

h 

— 1 

Z2 

X 

X 

Cv! 

X 

o 

y 

L_ 

• 

X 

O 

> 

{>' 

U 

X 

•*I 

• 

• 

• 

• 

• 

• 

I 

»-♦ 

<1 

t_ 

O 

U 

< 

c» 

C. 

r 

ri 

y 

Cj 

CD 

CT> 

cz 

L_ 

O 

1— • 

C 

.$ 

-J 

II 

c * 

_ j 

II 

'i- 

— S 

n 

r~ 

II 

ii 

11 

• 

ii 

II 

II 

_j 

K» 

ii 

II 

II 

II 

II 

II 

«i 

< > 

”■> 

v^' 

.» 

•s 

X 

2 

■S 

r-s 

Li. 

cn 

“5 

U_ 

IU 

t-4 

iH 

1-< 

»— v 

Q 

u- 

C/ 

pn 

•s 

X 

Uj 

•s 

“5 

2 

2 

X 

II 

II 

iH 

CV 

t') 

i-« 

(V 

O') 

2 

2 

CD 

X» 

•v 

»v 

2. 

X 

o 

M 

rv 

•v 

2 

rfT. 

i~4 

M 

~5 

IS 

»v 

•» 

rv 

rs 

rv 

II 

2. 

H- 

T-1 

r . 

X. 

w.. 

f» 

w 

h-l 

■~1 

r— 

r-H 

A- 

•“> 

~J 

“j 

J 

*“.> 

U 

rH 

2. 

•— 

w 

i — 

▼ 4 

2 

X 

I- 

1-1 

•s 

•v 

r- 

rv 

rv 

r> 

h 

Cr 

d 

2 

O 

«-i 

X 

O 

2 

V-1 

X* 

l-i 

«— t 

r— 1 

Xi 

»-4 

-J 

o 

2 

X 

o 

l_ 

O 

— 1 

5 

X 

G 

O 

C 

c 

W 

w 

C 

d 

a 

0: 

l: 

IH 

c 

c 

c 

CL 

c 

X 

X 

r* 

rv 

u: 

O 

bJ 

“•> 


cr  c 

( 

C 

*-<  O' 

<_ 


o 

c C 

4 < 

c 

C'  c 
t c 


O w>  o 
' c c. 
< ( ( 
t c_.  c 
('  r C' 
t c c 


o 

c 


90 


■ 


^ ■ - 


<Y(I«-1,J,‘43W>  - Y(I,J,NOW)  t /3Y 


oo  uj  co  oC)  cr>  o*>  cfj  a a to  co  o cr* 

K ^ k k K r K.  K.  ^ N.  K.  K 

S < * v * * * 

n a a x r 
XXXXXXXXT'TTZXX 
CO  O Cr  * 0 c 0 i 0 l ^ IV  iv  IV  IV  LC  CO 


HI  U\  >- 
w • o 

^ TH  Vy 

X I ^ 


^ O -r-» 
*5  *>  • 


IT 

3 

•k 

•k 

CJ 

HI 

O 

▼4 

rH 

w» 

w 

▼4 

X 

• 

1 

^-k 

V 

a 

4 

•k 

H4 

HI 

“0 

^■k 

n 

•> 

H 

X 

0 

4 

>- 

f^J 

▼H 

H 

H 

y-4 

—* 

4 

X 

f N 

4 

•» 

1 

1 

^k 

H^ 

>- 

I-; 

V— 

> 

4 

^k 

w> 

*■“* 

V) 

X 

4 

r.' 

HI 

• 

^ 

3 

a 

3 

o 

< 

<r 

V 

X 

O 

V 

O 

\ 

-J 

_i 

rvj 

H- 

rv 

a. 

1 

X 

X 

•k 

X 

n 

H* 

rr* 

K 

3 

c 

*. 

^k 

#k 

cr 

«T 

X 

W 

cr 

■V 

3 

^ t 

3 

•» 

4 

w 

•ft 

*. 

v_^ 

" r 

»k 

r-» 

4 

t— . 

«■ 

w 

* 

<- 

C. 

T-l 

7 

•»« 

7^ 

*k- 

*~  ■ 

^-, 

tH 

3 

X 

4 

»> 

>- 

1 

e 

HC 

•» 

^k. 

*) 

H 

4 

4 

4 

1 

H4 

mm) 

7-k 

▼4 

~N 

•k 

#. 

1 

> 

► 4 

^*» 

— > 

*k 

»■ 

— ) 

1 

»» 

H 

-o 

•* 

V 

•k 

> 

^-. 

«-i 

•k 

“N 

I-: 

4 

c 

«rs 

•* 

4 

^~s 

H4 

r 

H 

3 

1 

4 

«k 

M 

4 

M 

>- 

w 

4 

o 

•J 

H! 

f-4 

1 

HI 

X 

'W' 

• 

w- 

4 

#>« 

o 

> 

4 

• 

K, 

— 7- 

JL- 

+m* 

w 

1-4 

ZtZ 

> 

*•• 

h 

K 

j 

*. 

.3 

*k 

7 

7" 

w 

7- 

w 

O 

• 

1 

1 

7 

7' 

7 

4 

\ 

1 

fv 

3 

•k 

1 

T 

4 

•n 

w 

O 

H • 

w 

4 

*-k 

▼4 

V 

3 

7T 

3 

t 

4 

4 

3 

3 

3 

3 

4 

O 

•v 

3 

O 

V 

4 

4 

O 

<-> 

C 

k'l 

V 

“O 

(~5 

V 

X 

3 

7 

ZT 

3 

p 

•» 

7* 

•k 

• 

y-k 

■H 

*k 

•s 

•»  h 

tH 

o 

— > 

▼4 

*» 

CVJ 

HI 

*^l 

— N 

4 

tH 

•k 

-> 

r 

4 

—r 

4 

t4 

*k 

•* 

•k 

1 

v • 

1 

. 

4-1 

»« 

H 

HI 

4 

» -4 

II 

c 

M 

r-4 

»k 

‘H 

•k 

-> 

4 

”~ 

• 

W 

w 

7 

W 

H 

rH 

•k 

"N 

4 

•> 

y 

HI 

H 

1 

V- 

#k 

>- 

H 

tv 

y 

w 

7* 

w 

4* 

•k 

|H 

h.' 

w 

w 

w 

w 

^ 4 

>— 

H 

\ 

«T 

X 

7* 

X 

7" 

w 

*-< 

’ 

■v 

f i 

a 

i- 

0 

» 

1 

V 

1- 

_ 1 

*/ 

* 

> 

4 

K 

r' 

N 

N 

f' 

H 

ll 

H-l 

It 

II 

1-4 

II 

taJ 

o 

w 

HI 

w 

<a 

c 

CJ 

c. 

w 

»— - 

k_^ 

3 

>- 

rvi 

•k 

H- 

w 

>—* 

*-* 

*-H  tH  (\j 

4-  w X 

a hi  v x -o 

h w •> 

O X ^ — ‘ H< 

a 3 4 ir.  ^ 

4-  | C • X 

X ^ Cu  ^ 

4 »•  I | — .j 

HI  "7  •*  H4  (/i 


m 4 3 y 

^ *^  4 x r* 


-j  * 

4 4-  w 

^ ^ V 
* X ^ ^ 
.J.  c — ~ 

• X — > 3 

— • •»  •»  o 


^ > X rH 

— » w «-*  I 
4 w -J  Hi 
k-J  * ^ 

w ^ V 

O > 4 
X C 4 

tV  V -V 
Nk  ^ 

~ ~ 3 

O IH  C' 


H 4 _j  •» 
O HI  HI 


• M 

(\  ^ 

o - 

it  X n 
u w 

V,  v 

H-  H <T 
V C _J 
«i-  * 

X h 
► o -* 


ll  ll  it  ll  ► ♦"  . r-»  r— t 

w | W w I w 

r ^ r r a u a a_ 

r n.  h n h m fvj  c 

r o Cj 


14  14  1 


th  f\j  ro  -i  un  *x.* 


4 X «-  > 

S'  .i  w 


J + l)  t-  L A M ( T » J ) )»(A7P(I)  *•  n7°(T  ) ) )/(2.*D7)  SHNV3 


T 

I t'  CJ  f * C'  o 

C C'  c.  o c 

d 

a 


9li 


?i  1 = 1,  IOC, 


OUT? JT  1F?~M  7(0*!)  = ",  E12.7,///)  T? 


OJ 

tti 

a) 

CO 

CCi 

K. 

X 

K. 

K* 

K- 

K. 

£> 

n 

a 

a 

C 

CJ 

7 

ZT 

2 

a 

cr 

(Y 

IV 

(V 

CV 

Si  J 

ll 

-J 

_J 

<i. 

<1 

c_ 

o 

CO 

CO 

V 

\ 

3: 

O 

o 

U> 

2 

1- 

^ cC 

CV 

UJ 

Lu  (V 

L— 

l~ 

o 

l-« 

►-4  C* 

L- 

O' 

V »-- 

*? 

•»H 

T-l 

o 

— » 

CV 

“3  o 

CV 

CO 

rl 

II 

•*  c* 

r— < 

II 

_J 

fl 

> 

-J 

II 

o 

> 

w 

N1 

“3 

>- 

Osi  co 

• 

/V 

2 

• 

c* 

• 

*-4 

V 

• r 

r~> 

Li 

CP 

▼H 

OJ 

r 

cz  l. 

CV 

• 

tH 

• 

n 

H • 

• 

• 

h 

9 

h- 

1- 

r-s 

1 1- 

1- 

K 

O 

v*  1 

c 

> 

i-  c; 

1 

U> 

-J 

• 

• 

• 

"0 

• -j 

• 

• 

n 

^-1 

> 

> 

#• 

r~'  i_. 

N, 

r 

Ci' 

H 

~> 

I--4 

0 

—> 

II 

w 

W 

II 

W 

u 

> 

a 

Ll 

u 

U 

h 

h-< 

“0 

hH 

h~i 

*•* 

b- 

'—i  • H CVl  K;  O vx ‘ r*-  Ct  vJT  rH 

rr:  ::::::::: 

: ii  i*  ii  ii  n ii  ii  n ii  it 

T~«  II 

U ll  ^-1  H <»H  ^r-i  y~*  *~\  tH 

II'  V ▼ -»  *~i  *-»  r-i  r-"  *-*  r-1  ^r^  t-< 

2 <H  *«>*•»  * * ..  «,*, 

I 1 *—  ^ > s*  • l v»»  v.  w *W  * 0 ■ • . 

I-  CV  a H n M r*  ) . . u''  u f-  *.:•  .7* 

O O I!  —I  i f » I _J  „J  _J  _J  _j  _J 

C.  C‘  K *-»  C *.  ct  tj  c C tl  a tl 


9 1 NH 


T 


aj  co  ao 

K.  X X. 


a'j  m on  an 

k k n.  k. 


Q Q 0 O.  O ti  7 

7*  z?  7i  7 r 7 I 

cv  rv'  rv  <v  cv  rv'  to 


a>  cn  -O 

K.  k. 

X C X 
I7X 

10  (V  ir, 


r > 

r > 

«-4 

CD 

n 

CM 

*T 

_J 

e. 

CO 

d 

O 

h~ 

> 

CM 

Li 

z 

C 

O 

o 

X 

*"> 

X 

X 

2: 

s-N 

ci 

O 

o 

Or.' 

cm 

X 

CL' 

h-l 

C 

x. 

»~c 

u: 

*- 

X 

> 

•n 

<1 

o 

K 

C > 

IT. 

II 

X 

1 

o 

ll* 

1-4 

o 

I— 

_J 

H 

■f 

; 

: »> 

~y 

-J 

c~ 

cr  ^ 

ti  1 

a 

X 

2 

: -o 

i 

1 

X 

X 

O 

•* 

h- 

h~4 

O 

II 

n — i ^ 

♦— 

_J 

X 

<3 

00 

iLi 

to 

CZ- 

co 

«~* 

_J  C- 

< 

> 

1— t 

> 

H 

cm  d o 

_J 

X 

V »vM  ^ rl 

— V V * ~ 

ii  n n h cv  y 

1 I II  O H li 
C'  (/,  “o  T-i  17 


x ►-  i c 

r»  x r-. 

x u . jo  x i«j 

iT  4-  d ID 


o 

o 

CO 

h~ 

77. 

II 

II 

i:  z. 

a. 

7 

o 

1 

«j 

t~ 

li 

b- 

b r •. 

Cl 

CO 

-J 

_J 

CM 

i: 

X 

h- 

1- 

L- 

w ► - 

X 

U 

m 

ii. 

71 

•'O 

X X 

II 

ii 

C) 

C 

_J 

—1 

O 

o 

O 

tr 

C‘  li 

Li  O 

77 

Ui 

-J 

d 

c 

a 

u_ 

C-> 

— J 

A-  4. 

h 

en  (_, 

L 

CO 

O C"  c c 


X 


>5 

* 

CM 

O 

• 

4 

-J 

to 

O 

o- 

CO 

-J 

CD 

X* 

to 

L» 

o 

X 

tv 

XI 

o 

X— ' 

X 

T-f 

w 

4- 

w 

* 

L_ 

^x 

,-x 

*r 

O 

LD 

ro 

o 

Z' 

r*. 

ZF 

• 

• 

b 

<r 

to 

nl 

X 

cr* 

l > 

rv 

• 

:V 

rH 

<r 

X 

< 

o 

• 

4- 

_j 

_J 

.»• 

-J 

2_ 

J- 

4- 

n. 

00 

3 

^x 

rH 

M. 

• 

* 

cr. 

X 

^■x 

• 

!f\ 

CD 

CO 

H 

4- 

O 

CD 

i— < 

a , 

CD 

j 

m 

Ll. 

III 

O 

▼H 

CJ 

c_ 

Uk 

CD 

li\ 

o 

• 

(V 

• 

2L 

O' 

CM 

• 

■H 

y* 

o 

*r 

CM  -s 

O 

•“T 

• 

1TV 

1 

1 

o 

CM 

O 

X 

CO 

o 

CM 

tu 

U- 

CM 

x-» 

4-  •*- 

u 

X 

fo 

K. 

X 

V 

1 

li 

<o 

| 

X 

V 

III 

*— i 

lx 

rr. 

r> 

s* 

IA 

• 

V- 

«A 

n 

i*' 

3“ 

CD 

x- 

•r 

#"X 

L 

CJ 

w 

• 

1 

if 

X 

O 

4 

* 

CM 

c 

V 

a 

CM 

r 

•f 

fl 

O 

V 

UJ 

* 

cv. 

w* 

o 

o 

C ) 

£~ 

V 

• 

10 

*» 

D 

r 

X 

-J 

V 

o 

CM 

10 

X 

IU 

UJ 

CD 

CJ 

<r 

rv 

co 

1-4 

Z 

o 

or 

& 

0- 

X 

lx 

> 

<3 

Cv 

K 

cr 

o 

> 

Q' 

CL 

ll 

V- 

<1 

V 

X 

C 

IT 

If* 

21 

•— 

•X 

_l 

_J 

4- 

rv 

CP 

00 

0 

c 

h-i 

1 

to 

cv 

Xt 

V- 

V 

w- 

L' 

•— 1 

.. 

V 

CM 

V 

Cl 

• 

i:\ 

ar ; 

X_*. 

N. 

c 

V 

<_ 

cH 

ll 

• 

• 

X-* 

r' 

e 

rr 

o 

• 

Zt~ 

x^ 

r» 

o 

_J 

U! 

<z 

►fc 

*-* 

• 

n. 

1/ • 

ro 

to 

CJ 

c 

II 

• 

1 

• 

O' 

/-X 

II 

(ii 

4 

X 

4- 

rr> 

ry 

V 

h- 

• 

< 

• 

( : 

l 

cm 

X 

^■x 

o 

C 

i'- 

l: 

C 

Q' 

a' 

1* 

Ll 

fD 

UJ 

3 

3' 

L 

a: 

If 

\iD 

•'C 

^ 

• 

ll 

2 ' 

iA> 

<^1 

<-* 

Cv. 

•w 

f*. 

> 

H 

> 

h- 

X 

Cr' 

• 

X 

• 

> 

Cxf 

K 

>- 

CV' 

< 

_l 

x-' 

f 

21 

x-' 

*0 

it 

0' 

* 

ll 

<x 

O 

0 

X 

1- 

c 

hH 

V>  ('• 

#A 

O c • *4  *4 

^ • • • 

X l*.  llJ 

> c c e 

• • • 

li.  li  U 
II  X X * 
> > > 
lJL  w w w 

> U li  1 

> H I i X 


C t.'  c c 


C‘  c c 


u . 

c: 

ir\ 

CD 

|> 

cv 

X 

r* 

• 

1 

N 

• 

i 

• 

• 

L 

c 

cr 

cr 

O 

CM 

u. 

CD 

CD 

x 

f- 

cr 

1 r 

ir 

V 

*-» 

X 

«v 

<r 

^-x 

w» 

x^ 

w 

ir 

X 

X 

X 

cv. 

^-x 

«-x 

*-X 

• 

• 

• 

• 

• 

it 

%'x. 

O' 

1*0 

t* 

w 

y* 

1 

r.j 

1 

1 

1 

1 1 * 

? _ 

• 

• 

• 

• 

cr 

4- 

r 

CO 

CM 

• 

«r i 

tH 

U ' 

w 

w 

x_ 

o 

• 

II 

II 

ro 

II 

x 

II 

ll. 

X 

Ia. 

M 

U- 

e 

c- 

c 

AA. 

X 

y 

H 

X 

y 

X 

V 

x^ 

xj 

< 

t 

«"4. 

IL 

lu 

C 

ll 

cr 

j 

o 

li 

f 4 

1‘ 

or 

r. 

C 

c. 

a 

c 

c. 

C- 

1 

i 

o 

f 

X 

CC 

CC  cr  K 


rins  I?  LIME**  CTT  TO  :~LbSSrl')~  FOR  =>A  5TICJL  A9  YIELO 


Appendix  G 


Estimate  of  Apparent  Airshock  Velocity 
Along  Ground  Surface 

The  shock  front  from  a nuclear  hurst  can  be  considered  a 
spherical  shell  which  is  expanding  with  some  average  velocity, 

V . The  intersection  of  this  spherical  shock  front  with 

ci V G 

the  ground  surface  forms  a circular  region.  This  circular 
line  of  intersection  will  expand  along  the  ground  with  a 
velocity,  called  the  apparent  velocity  V , which  will  differ 
from  Vavg.  If  the  airshock  is  approximated  locally  as  a planar 
wavefront,  the  following  figure  describes  the  position  of  the 
airshock  at  two  times,  t and  t+dt. 


In  dt,  the  shock  wave  radius  S (also  equal  to  the  slant 
range)  expands  by  Vavedt.  Thus, 


S( t+dt)  = S(t)  + V dt. 

ave 


Also , 


(dr)^  = (V^  dt )‘ 

ave 


S(t)  d(f)} 


, or 


— = V^  + 

dt  ave 


S(t) 


1 


dt 
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The  surface  velocities  resulting  from  simulation  of  a 28 
kiloton  atmospheric  hurst  at  500  feet  height  of  burst  yielded  a 
dust  layer  ballistically  reaching  only  0.64  cm  at  its  highest 
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induced  precursor  as  postulated  (ballistic  rise  only)  fails  to 
re-create  the  2-3  meter  high  dust  layers  observed  in  experi- 
mental atmospheric  nuclear  testing.  However,  the  motions  are 
felt  to  be  significant  enough  to  be  included  in  any  attempt  to 
model  from  first  principles  the  precursor  and  the  up-sweep  of 
dust  behind  the  shock  front. 


